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Abstract 

A developable surface can be formed by bending or rolling a planar surface without stretch- 
ing or tearing; in other words, it can be developed or unrolled isometrically onto a plane. 
Developable surfaces are widely used in manufacturing with materials that are not amenable 
to stretching. A ship hull design entirely composed of developable surfaces would greatly 
reduce production costs of that hull. 

This thesis describes a new, user-friendly method of designing developable surfaces with 
a B-spline representation. First, it expands the work of Aumann in designing developable 
B-spline strip surfaces whose directrices lie on parallel planes. The computer program 
developed can assist in designing surfaces from curves of any degree with any number of 
segments, and includes a test for regularity. Second, a new method is developed which 
permits the design of developable surfaces with general three-dimensional space curves as 
directrices. The computer program permits design of degree (3-1) B-spline developable 
surfaces with two patches using a minimization process. The basis is provided to extend 
this to surfaces with more patches and higher degrees. A test is provided to ensure that the 
minimization process results in a developable or nearly developable surface. This method 
is extended to include special interesting cases such as triangular degenerate patches and 
surfaces with a planar resulting directrix. 

This thesis also treats common differential geometry properties such as lines of curvature 
and geodesics that are useful in the design and manufacturing process. Lines of curvature 
are vital in the forming process, since the planar shape must be placed so the rollers are 
parallel to the of zero curvature. As an inflection line greatly affects the forming 
process, a method is aescribed to determine the inflection line in advance; several properties 
of developable surfaces related to inflection lines are also described. 

Straight lines on a plane map to geodesics on a developable surface. This fact is of 
assistance in many aspects of the use of developable surfaces, including layout and quality 
control. A method is described to calculate geodesics on a developable surface as an initial 
value problem instead of the more complicated boundary value problem. 

Engineering examples are provided for each topic, and a small boat hull is designed 
using the methods described. In addition, recommendations for further research are made. 
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Chapter 1 

Introduction 



1.1 Background and Motivation 

A ruled surface is a curved surface which can be generated by the continuous motion of 
a straight line in space along a space curve called a directrix. This straight line is called 
a generator, or ruling, of the surface. Developable surfaces are a subset of ruled surfaces 
which have the same tangent plane at all points along the generator. A developable surface 
can be formed by bending or rolling a planar surface without stretching or tearing; in other 
words, it can be developed or unrolled isometrically onto a plane. Developable surfaces are 
also known as singly curved surfaces, since one of their principal curvatures is zero. 

Developable surfaces axe widely used in manufacturing with materials that are not 
amenable to stretching. Applications include the formation of ship hulls, ducts, shoes, 
clothing and automobile parts such as upholstery and body panels [14]. 

In shipbuilding, developable surfaces are shaped using only rollers or presses. Heat 
treatment is then used only to remove distortion induced by welding or other means. Dou- 
bly curved surfaces, on the other hand, must be heat treated after rolling to induce the 
additional curvature. The heat treatment is normally done by hand, by a skilled artisan 
with years of training to achieve the correct amount of bending. This is an extremely time 
consuming, labor intensive and thus expensive process. 

According to Avondale/IHI Shipbuilding Technology Transfer data for a tanker, only 
15.1% of the curved plates in a ship hull are singly curved, while 65.8% of the plates are 
doubly curved, requiring roller and heat treatment processes [28]. An intensive effort to 
increase the amount of developable surfaces in the hulls of merchant ships at Burmeister 
& Wain Shipyard [34] has resulted in a reported 20% reduction in manhours required to 
produce a hull. Designing a ship entirely of singly curved, or developable, surfaces would 
reduce construction costs even more. 

Recently researchers in Computer Aided Geometric Design have been quite active in 
investigating the representation of developable surfaces in terms of Non-Uniform Rational 
B-Spline (NURBS) surfaces or a special case of NURBS called Bezier surfaces [1, 3, 14, 
29]. NURBS curves and surface patches are the most popular representation method in 
CAD/CAM due to their generality, excellent properties and incorporation in international 
standards such as IGES (Initial Graphics Exchange Specification) and STEP (Standard for 
the Exchange of Product Model Data). Thus, it would be beneficial to design developable 
surfaces using a B-spline representation. 
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1.2 Research Objective 



Recently, efforts have been underway to revitalize commercial shipbuilding in the United 
States. Several professors at MIT are involved with developing cost saving methods in ship 
fabrication to ensure that shipbuilding in the US will be competitive in the world market. 
This research is in support of one of those efforts. 

The main goal of this research is to develop a user-friendly method of designing devel- 
opable surfaces with a B-spline representation. The effort is then extended to address some 
common differential geometry properties that will be useful in the design and manufacturing 
process. 

The ultimate goal is to provide a method to design a complete ship hull from developable 
surfaces and to generate cutting and bending information in a format that is user friendly 
for both the engineer and the worker. Although this thesis does not go that far, it takes a 
major step toward this goal. 

1.3 Thesis Organization 

In this thesis, Chapters 2 through 4 review basic differential geometry properties and intro- 
duce the concepts of developable surfaces. They also briefly review the Bezier and B-spline 
representations of curves and surfaces. The information introduced in these chapters will 
be used throughout the thesis. 

Chapter 5 describes a new user friendly method for design of developable surfaces in 
B-spline representation. Section 5.1 reviews the current literature on developable surfaces. 
Section 5.2 describes the design of strip surfaces that are constrained between two parallel 
planes. Section 5.3 describes the design of surfaces with directrices that are space curves 
and treats special cases such as triangular degenerate patches and surfaces with a planar 
resultant curve. Section 5.4 describes the unrolling, or development, of the surfaces into a 
plane. Section 5.5 gives some engineering examples. 

In Chapter 6, a geodesic on a developable surface is found as the solution to an initial 
value problem rather than a two point boundary value problem. The geodesics can be used 
in ship design for laying out butts and strakes, among other uses. In Chapter 7, the lines of 
curvature on a developable surface are analyzed and the line of inflection is deflned. These 
concepts are required for the bending of steel plates into developable surfaces; the steel 
must enter the rolls in a direction parallel to the lines of zero curvature, and cannot be 
rolled past a line of inflection [33]. 

A small boat was designed using the methods described in this thesis. The boat is 
presented in Chapter 8 as a practical engineering example. 

The thesis concludes with Chapter 9 which includes recommendations for further re- 
search. 
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Chapter 2 

Review of Curves and Surfaces 



2.1 Introduction 

This chapter reviews the basic theory of curves and surfaces including such topics as the 
Serret-Frenet formulae, the first and second fundamental forms, curvatures and geodesic 
curves which will be used in derivations of many formulae in the later chapters. The 
information in this section follows the derivations in classical texts on differential geometry 
and geometric modeling such as those by Kreyszig [26], Struik [41] and Faux and Pratt [11]. 



2.2 Basic Theory of Curves 



Throughout this thesis, curves are represented in a parametric manner, r(u), as a function 
of one parameter, u, that lies within a closed interval u\ <u<U 2 - The curve is a mapping 
from a one parameter interval to three-dimensional space 

r(u) = [i(u) y{u) z{u)]'^. 



The points on a curve are regular as long as at least one of the first derivatives is not equal 
to zero. In other words, a point is singular if dx/du = dy/du = dz/du = 0. The curve can 
be reparameterized if u can be expressed as u = f{u\) as long as du/du\ ^ 0. 

Arc length is defined ais the distance along a curve between parameter values uq and u, 
and can be represented as 



(u) = f yjx^ + ip’ + du = f 

JuQ JuQ 



(2.1) 



In this thesis, derivatives with respect to the arc length s will be represented by a prime 
and derivatives with respect to u will be represented by a dot. Derivatives of arc length s 
with respect to parameter u and vice versa are 



s 

s 



ds n — r 

— = vr • r = 
du 

ds r • f 
du y/r t 
du 1 

ds |r| 
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du' _ r r 
ds (r • r)^ 



The unit tangent vector t to a curve is defined by the unit vector that passes. through 
two points in the curve, r(s) and r(s + h) (see Figure 2-1), as h approaches zero; i.e. 



t 



r(s -I- /i) - r(s) 
lim ; 

h — ^0 #l 

dr{u) du r{u) 

du ds lr(ii)| 



dr 

ds 



= r'(5) 



This final representation clearly shows that t is a unit vector. The line passing through the 
point r(n„) in the direction of t(nn) is the tangent to the curve at 




Figure 2-1: Tangent Vector 



All vectors normal to t{un) at Un lie in a plane called the normal plane. Much as 
the tangent is defined by two points that approach one another on a curve, the osculating 
plane is defined by three points that approach one another. In other words, the osculating 
plane is the plane in which the curve lies at a point on the curve. If the curve is planar, 
the entire curve lies in a single osculating plane, and the osculating plane is constant for 
the entire curve. To mathematically define the osculating plane following the derivation in 
Kreyszig [26, p.31], take three points at parameter values n, u -h /ii and u + /i 2 . The plane 
can be defined by the two linearly independent vectors a{ = r{u + hi) — r(w), i = 1,2 or 
linear combinations thereof. Let 

(i) ^ r(u -t- hi) - r(u) 
h^ 



and 



2(y(2)-v(l)) 

/i2 — hi 
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The Taylor series expansion of r{u + h^) is 



/j2 

r(u + hi) = r{u) + hiir{u) + -^r(u) + o(/i?) 



( 2 . 2 ) 



where o(h") is a vector of Landau symbols o(h”) with the property that [41, p.3] 

lin,2^ = 0. 

h-^0 /l" 



Using (2.2), we obtain 

= 



and 



r(u + /ii) — r(u) 

r(u) + hiir{u) 4- ^hir{u) + o{hi) - r{u) 

Ih 

./ . 1 , \ o(^i) 

r(u) + -hir{u) + 



w = 



2 (v( 2 ) - v(i)) 



h2 ^ h\ 
2 



fi2 — h\ 



r{u) + ^h 2 r{u) + 



- fr(«) + \hir{u) + 



= r{u) 



2{hio{h2) — h2o{h\)) 
hih2{h2 — hi) 



As hi 0, r(u) and w — r(w). Therefore, the osculating plane is defined by r(i/) 

and r{u). Note that r{u) must be at least twice continuously differentiable. 




Figure 2-2: Osculating Plane 

The principal normal is defined as the intersection of the osculating and normal planes. 
Using an arc length parameterization, we find that differentiating t • t = 1 with respect to 
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s yields 



t • t' = 0, 



so t' is perpendicular to t and lies in the normal plane. Since t = r' = ru', 

t' = ru'^ + r(u')^. 



Thus t' lies in the plane of r and r, and therefore in the osculating plane. Since t' lies in 
both the osculating and normal planes, a unit vector in the direction of t' is the unit normal 
vector 



t' 




The magnitude of t' is called the curvature and is represented ets ac. The curvature vector 
k is represented as 



k = t' = Kn. 



(2.3) 



A physical sense of the curvature can be found from the simple example of a line, which 
has a constant tangent vector, so t' is zero. Therefore, « = |t'| = 0, and a line has zero 
curvature. Parameterized by u, k and k can be represented as [41, p.l6] 



k = 



K = 



A.f •//,••/ f\2 

t = ru r[u ) = 
(Acn) • («n) = 



(r • f)r — (r • f)r 
(r • r)2 



‘(r • r)f — (r • f)r' 




'(r • r)r — (r • f)f ' 


(r X r) • (r X r) 


(r-r)^ 




— 1 
<N 


(r • r)^ 



where the Lagrange identity (a x b) • (a x b) = (a • a)(b • b) — (a • b)2 is used. 
The normal to the osculating plane is the binormal, which is defined as 



b = t X n. 



Since these three vectors are unit vectors that are mutually orthogonal, they satisfy the 
relations 



bb=l nn=l tt=l 
b-n = 0 n-t = 0 t-b = 0. 

Differentiating b • b = 1 and b • t = 0 yields 

b' ♦ b = 0 



and 



b' • t = — b • t' = — b • Acn = — Ac(b • n) = 0. 

Therefore, b' is orthogonal to both b and t and must then be parallel to n. Let 

b' = — rn 
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where r is the torsion of the curve, A physical sense of torsion can be found from the fact 
that a curve with zero torsion and nonzero curvature is a planar curve. 

Finally, an equation for the derivative of the normal vector corresponding to the rate of 
change of the tangent plane [26, p.41] can be determined. Differentiating n • n = 1 yields 
n • n' = 0. Therefore, n' is orthogonal to n and must satisfy the relationship 

n' = at+/?b. (2.4) 

Multiplying equation (2.4) by t and b respectively yields 

a = n' • t and ^ = n' • b. 

Differentiating n • t = 0 yields 



a = n' • t = — n • t' = — n • nn = —k. 
Similarly, differentiating n • b = 0 yields 

/? = n' • b = — n • b' = — n • (— rn) = r. 



Thus, 



n' = —Kt + rb. 



The representations of the derivatives of the three defining vectors for a curve are termed 
the Serret- Prenet formulae [41, p.l8] 

t' = «n 
b' = — rn 

n' = —Kt + rb, 

which describe all aspects of a space curve, as shown in Figure 2-3. 




Figure 2-3: Space Curve 
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2.3 Basic Theory of Surfaces 



A general parametric surface can be defined as a vector-valued mapping from a two- 
dimensional parametric ui;-space to three-dimensional space 

R(u, u) = [x(u, u) y{u,v) z(u,u)]^, 



where the two parameters u and v lie in the closed intervals u\ <u <U 2 and v\ < u < V2- 
The surface is regular at all points where the matrix 



dx 


dy 


dz 


du 


du 


du 


dx 


dy 


dz 


dv 


dv 





(2.5) 



has rank 2, or equivalently x Ru 7*^ 0 [41, p.56]. Partial derivatives of a surface will be 
represented as 



^ 3R , „ 3R 
= — and R^ = — . 
ou ov 



A curve on the surface is represented as r(u, v) with u = u(^) and v = v{t) or equivalently 
R(u(t), u(^)) = r(t). This occurs when the rank of the matrix (2.5) is 1 at every point. 
When we keep u constant in R(u,u) by setting u = Un, we have an equation for a curve 
that depends only on u, and is thus called an isoparameter curve. Similarly, u = is also 
an isoparameter curve. 

The vector 



dR _ du dv 

dt dt dt 



(2.6) 



is tangent to the curve when u and v are functions of t. This is of course also a tangent to 
the surface. The tangent plane to the surface at any point is determined by the two vectors 
R^ and R^. The normal to the surface is orthogonal to this plane. A unit vector in the 
direction of the surface normal is given by 



R^i X R^ 
|Ru X R^ 



(2.7) 



Note that this definition of the normal vector is undefined when R^ x R^ = 0, supporting 
the assertion that a point on the surface is singular if R^ x R^ = 0. This may be due to 
the shape of the surface or the choice of parameterization. In Section 5.3.5, an alternative 
definition of N is explored for degenerate points. 

Using the definition of arc length (2.1), the distance between two points on a curve on 
the surface is found by integrating 

ds = y/ dxdx -h dydy -h dzdz = VdR • dR. (2.8) 



Combining (2.6) and (2.8) yields the first fundamental form of the surface 

/ = ds^ = {Ru,du -h Rvdt^) • {Rudu + Rydv) = Edu^ -{- 2Fdudv -h Gdv^ 
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where 



Note that 



E — • Rix, 

EG-F^ = 



F = Ru ■ RtM G = R,; • R,;. 



(Ru • R^)(R^ • Rv) — (R^ • Rd)^ 
(Ru X R„) • (R^ X R„) 
|(R„xR„)|2. 



Thus, EG - F2 is always positive if the surface is regular. 



(2.9) 



The first fundamental form enables us to determine arc lengths, angles and areas on the 
surface. In order to determine curvatures, the second fundamental form is required. First, 
the curvature vector of the curve on the surface, found from equation (2.3), is split into its 
normal and tangential components 

dt 

k = — k,2 “ 1 “ k^ — “ 1 “ kp 

where k^ is the normal curvature vector and Kn is the normal curvature of the curve on 
the surface at that point. The tangential portion of the curvature vector, k^, is called the 
geodesic curvature vector. It will be discussed in more detail later. 



If the angle between the surface normal and the normal to the curve is designated 7 as 
shown in Figure 2-4, then [26, p.118] 

n • N = cos( 7 )|n||N| = cos( 7 ) 

from the definition of a dot product. Since t' = R" = k • n, 

KCOs{j) = K • n • N 

= R"-N 

= (Rudu -f Rydvy • N 

= (Ruududu -f 2Ruvdudv -f Ryydvdv -f Ryd^u -f Ryd^v) • N. (2.10) 

Since R^ . N = R^; • N = 0, equation (2.10) reduces to the second fundamental form, 
represented by 

II = {Ruududu 4- 2Ruydudv + Ryydvdv) • N 
= Ldudu -f 2Mdudv + Ndvdv 

where 



= Ruu • N 
= Rut; • N and 
= Ryy • N. 



L 

M 

N 
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( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 



Also note that since • N = 0, 



and therefore 



Similarly, • N = — 



d(Ru • 
du 


= R.. 


•N + : 


Ru • Nu 


= 0, 






II 


-Ru 


•Nu. 






Ru • N„ = 


- * N^-y and 


•N = - 


-Rt, • Nt,. Therefore, 




L = 


-R„ • N„ 








(2.14) 


M = 


-Ru • N„ = 


-Rt,- 






(2.15) 


= 


-^(Ru-N„ 


+ Rt; 


•Nu) 


and 




N = 


-R„ • N„. 








(2.16) 



N 




Figure 2-4: Curvature Vectors 



To find the normal curvature, we begin with the fact that the surface normal is orthog- 
onal to the tangent to the curve. Differentiating t • N = 0 with respect to s yields 



dt dN 



Thus we have 



Therefore, 



dt 

k.N = — -N = 
ds 




dR dN 
ds ds * 



— k N — ^ ^ ~ (Rudu -1- Hydv) • (Nudii -1- Nydv) _ II 
ds ds ds‘^ I 

_ Ldudu + 2Mdudv H- Ndvdv 

Edudu + 2Fdudv -h Gdvdv 
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L + 2MA + NX^ 
E + 2FA + GA2 



(2.17) 



where A = dv/du. The negative sign in equation (2.17) ensures that if is positive, the 
center of curvature lies opposite to the direction of the surface normal. 

Any point on the surface has many curves passing through it, each of which has an 
associated The majcimum and minimum values of Kn at a point are called the principal 
curvatures of the surface at that point and are designated k\ and « 2 - The majcimum and 
minimum values of occur when ^ =0, [11, p. 112] 

dK 

^ {E + 2FX + G\^){2M + 2NX) - {L + 2MX + NX^){2F + 2GX) = 0. (2.18) 

dX 

Using (2.18) in (2.17) yields 

L + MA _ M + NX 
- E + FX~ F + GX 

given 

{E + 2FX + GX^) = {E + FX) + X{F + GX) and 
(L + 2MX + ATA^) = {L + MX) + X{M + NX). 

Therefore the extreme values of Kn satisfy the two simultaneous equations 

{L KnE)du {M KnF)dv = 0 and (2.19) 

(M + KnF)du + (AT + KnG)dv = 0. (2.20) 



These equations can be simultaneously satisfied if and only if 



L 4" KnF Af -j“ KnF 
M 4- KnF N KnG 



The discriminant of this quadratic equation in Kn is greater than or equal to zero. Therefore 
the equation has either two distinct roots Kmax a-nd Kmin, the maximum and minimum 
principal curvatures, or double roots Kn^ the normal curvature at an umbilical point. The 
corresponding directions A define directions in the uv-plane. The corresponding directions 
in the tangent plane are called principal directions of curvature and are always orthogonal 
except at umbilical points^ where Kmax and Kmin are identical. In the special case where the 
identical principal curvatures vanish, the surface becomes locally flat. The two roots are 
given by 



f^max = H - K and 

Kmin = H - - K 

where K is the Gaussian curvature and H is the mean curvature defined by 

_ LN-M^ 

- eg -F^ 

__ 2FM -EN-GL 
2{EG - F2) • 



( 2 . 21 ) 

( 2 . 22 ) 



(2.23) 

(2.24) 
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Using (2.21) and (2.22) in (2.23) and (2.24), it can easily be shown that 



K 

H 



f^max^min 
1. 

‘^max H“ ^min) 



2(«n 



Note that at an umbilical point, = K and at a flat point, K = H = 0. 

The normal curvature can be represented in terms of the principal curvatures k\ and «2 
■ as 



Kn = «i cos^(a) + K 2 sin^(a) (2.25) 

where a is the angle between the direction dv/du and the direction dv = 0. This is known 
as Euler theorem [26]. 

Recall that the curvature vector of a curve on the surface can be split into its normal 
and tangential components 



k — + kp, 

where the tangential component, k^, is known as the geodesic curvature vector. Let u be 
a unit vector in the tangent plane perpendicular to t such that t, u and N form a right- 
handed coordinate system in that order as shown in Figure 2-4. The geodesic curvature Kg 
can then be found from 



k^ = KgU. 

Additionally, since the magnitude of can be represented as |«cos 7 |, the magnitude of 
Kg can be represented as |Ksin 7 |. 

Geodesic curvature can be represented using only the first fundamental form. Since 
u • N = 0 and u • u = 1, 

= U • k = u • ^ = u • t' = (N X t) • t' = (tt'N). (2.26) 

ds 

The unit tangent vector of the curve on the surface is given by 

t = Ruw' + 

thus 



Therefore equation (2.26) can be rewritten as 
Kg = (t X t') • N 

= [(R„u' + R„t;') X (R„u(u')^ + 2R„„u't;' + Kvv{v'? + Ruu" + R„t;")] ' N 

= [(Ru X Ruu)(w')^ + (2Ru X Ru„ + Rt, X Ruu)(u')2u' 

+(Ru X R„„ + 2R„ X 'R^y)u' {v')"^ 

+(R„ X R„„)(t;')^ + (Ru X K,){u'v" - u"v')\ ■ N. (2.27) 



22 



Each of the coefficients of the combinations of u and v derivatives in (2.27) can be written 
in terms of the coefficients of the first fundamental form and their derivatives. For example, 



(Ru X Ru„) ■ N 



(Ru X Ru„) • (Ru X R„) 

VEG - F2 

(Ru • Ru)(R^u • Rt;) — (Ru • Ruu)(Ru • Rd) 
y/EG - F2 

2(£;(Ru„ • R„) - (R„ • B^u)F)VEG-F^ 

2{EG - F2) 

{2EF^ - EEu - FEu)y/EG-F^ 

2{EG - F2) 
rfiVFG-F2 



where Ffi is one of the Christoffel symbols defined by the coefficients of the first fundamental 
form E, F and G and their derivatives, F„, F„, G„, F„, F„ and G„ 

2 _ 2FF„ - FF„ + FEu 
“ “ 2(FG - F2) 

2 ^ EGu - FEy 

2(FG - F2) 

2 _ EGy — 2FFy + FGu 
22 “ 2(FG - F2) ■ 

Similarly reducing all the other coefficients in (2.27) yields 

= [r2i(u')' + (2r22-r}i)(u')2|;' + (r22-2r}2)u'(u')2-ri2(t;')' 

+u'v" - u"v'] VEG - F2. (2.28) 



GEy - 2FFu + FEy 

2{EG - F2) 

GEy - FGu 
2(FG -F2)’ 

2GFy - GGu + FGy 
2{EG - F2) 



A geodesic is defined as a curve with zero geodesic curvature [41]. Straight lines on a 
surface are geodesics since the curvature vector k vanishes. For geodesics that are curved, 
the curvature vector coincides with the surface normal vector and, since the curvature vector 
lies in the osculating plane, the osculating and tangent planes are normal. The equation for 
the geodesic can be obtained by setting Kg equal to zero in equation (2.28), which yields 

uV - u'v" = r 2 i(u ')3 + (2f 22 - rli)(u') 2 u' + (rl2 - 2r\^)u'{v')^ - ri2(t;')" 



assuming that the geodesic is everywhere regular so EG — is always positive. An al- 
ternative representation for a geodesic is given by the set of coupled second order ordinary 
differential equations [41] 



d^u 

d?v 

ds'^ 






11 



5i + r2^ 




du dv 

_9 du dv _9 




= 0 
= 0 



(2.29) 

(2.30) 



where the two equations are related by the condition ds^ = Edv? + 2Fdudv + Gdv^. 
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Chapter 3 

Curve and Surface Representation 



3.1 Introduction 

This chapter reviews the representation of curves and surfaces in Bezier and B-spline forms 
and treats the special properties associated with each. In addition, algorithms for knot 
insertion and removal, curve evaluation and splitting are summarized. The descriptions are 
based on texts by Hoschek and Lasser [22], Piegl and Tiller [36] and Yamaguchi [44]. 



3.2 Bezier Curves and Surfaces 

Bernstein Polynomials 

The Bernstein polynomials are defined as 

n! 



Bi,n{u) = 1 i)! ~ * = 0,...,n, 



(3.1) 



and have several properties of interest. The property of positivity states that each polynomial 
factor is non-negative such that 

Bi^nM >0, 0<u< 1 

for all i and n. The partition of unity property states that the Bernstein polynomials sum 
to 1 for all 0 < u < 1, or 

E5.,n(u) = l. 

1=0 

The derivative of a Bernstein polynomial is 

d>Bj^fi{u) 



du 



— ^[-Si— l,n— 1 (^) 



where = ^n,n-i = 0. The linear precision property of a Bernstein polynomial 



i=0 ” 
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allows the monomial u to be expressed as the weighted sum of Bernstein polynomials with 
coefficients evenly spaced in the interval [0,1]* 

Bezier Curves 

A Bezier curve is a spline curve that uses the Bernstein polynomials as a basis. A Bezier 
curve of degree n (order n -h 1) is represented by 

n 

*■(“) = XI 0 < u < 1. 

1=0 

The coefficients, bj, are the control points that determine the shape of the curve. Lines 
drawn between consecutive control points of the curve form the control polygon, A cubic 
Bezier curve is shown in Figure 3-1. Bezier curves have the following properties: 

• The first and Icist control points are the endpoints of the curve. In other words, 
bo = r(0) and b„ = r(l), 

• The curve is tangent to the control polygon at the endpoints. 

• The convex hull property states that the entire curve is contained within the convex 
hull of the control points. 

• The first derivative of a Bezier curve is represented by 

dr(u) 

r(u) = = « X (u), 0 < « < 1. 




Figure 3-1: A Cubic Bezier Curve with Control Polygon 

A Bezier curve can be evaluated at a specific parameter value no and the curve can be 
split at that value using the deCasteljau algorithm^ where [22, p.l25] 

bf(wo) = (1 - tio)bfr/ +uobf“\ k = l,2,...,n, i = k,...,n 
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is applied recursively to obtain the new control points. The algorithm is illustrated in 
Figure 3-2, and has the following properties: 

• The values are the original control points of the curve. 

• The value of the curve at parameter value no is b^. 

• The curve can be represented as two curves, with control points (bp , b^ , . . . , b^) and 

(K, b°). 

• Since this is merely a change in the basis representation, the curve itself remains 
parametrically and geometrically unchanged. 




Figure 3-2: The deCasteljau Algorithm 



Bezier Surfaces 

A tensor product surface is formed by moving a curve through space while allowing defor- 
mations in that curve. This can be thought of as allowing each control point to sweep a 
curve in space. If this surface is represented using Bernstein polynomials, a Bezier surface 
is formed, with the following formula: 

m n 

R(u,v) = ZE 

t=0 j=0 

Here, the set of lines drawn between consecutive control points b^j is referred to as the 
control net. An example of a bi-quadratic Bezier surface with its control net can be seen in 
Figure 3-3. The following conditions apply: 

• The boundary isoparametric curves (n = 0, u = l,u = 0 and ^; = 1) have the same 
control points as the corresponding boundary points on the net. 

• The corners of the surface coincide with the corner points of the net, and the deriva- 
tives in the u and v directions at the corners are tangent to the net. 

• Bezier surfaces exhibit the convex hull property. 
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Figure 3-3: A Bezier Surface with Control Net 



Continuity Conditions 

Bezier curves can represent complex curves by increasing the degree and thus the number 
of control points. Alternatively, complex curves can be represented using composite curves, 
which can be formed by joining several Bezier curves end to end. If this method is adopted, 
the continuity between consecutive curves must be addressed. 

One set of continuity conditions are the geometric continuity conditions, designated by 
the letter G with an integer exponent. Position continuity, or continuity, requires the 
endpoints of the two curves to coincide, 

r(i)(l) = r(2)(0). 

The superscripts denote the first and second curves. Tangent continuity, or continuity, 
requires the tangents of the curves to be in the same direction, 

r^^^(l) = ait 

r(2)(0) = Q2t 

where t is the common unit tangent vector. G^ continuity is important in minimizing stress 
concentrations and preventing flow separation. 

Curvature continuity, or G^ continuity, requires the center of curvature to move contin- 
uously pcLSt the connection point, 

r(2)(o) = r^^^(l) -h/ir^^^(l). 

G^ continuity is important for aesthetic reeisons and for preventing fluid flow separation. 
More stringent continuity conditions are the parametric continuity conditions, where 
continuity requires the kth derivative of each curve to be equal at the joining point. In 
other words. 



du^ du^ 
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The and continuity conditions for consecutive segments of a composite degree n 
Bezier curve can be stated as [44, p. 215] 



bni-l H 



(bfii bjii—i) — hi (byif-i-i bfii) ? and 

bni— 2) — + 7 — b7ii_^2) 



H-hl 



(3.2) 

(3.3) 



where, for the zth Bezier curve segment with u values between knot values ti and ti^i 
inclusively, hi — ^14.1 — U, See Figure 3-4. 




Figure 3-4: Continuity Conditions 



3,3 B-Spline Curves and Surfaces 

The Bezier representation has two main disadvantages. First, the number of control points 
is directly related to the degree. Therefore, to increase the complexity of the shape of the 
curve by adding control points requires increasing the degree of the curve or satisfying the 
continuity conditions between consecutive segments of a composite curve. Second, changing 
any one control point affects the entire curve or surface, making design of specific sections 
very difficult. These disadvantages are remedied with the introduction of the B-spline 
representation. 

A B-spline (basis-spline) curve is a piecewise polynomial curve meaning it is made up of 
polynomial curve segments. It is a spline curve with a different basis than the Bezier curves, 
although a Bezier curve is a special case of a B-spline curve. The B-spline representation 
has the advantage that the control points affect the curve only locally, and a B-spline of 
order k may have as many control points as required to describe the curve. 

A B-spline curve is a series of polynomial segments of degree A; — 1 joined together at 
knots tk- A B-spline of order k (degree A; — 1) with n + 1 control points has the equation 

n 

*■(«) = n>k-\, u e [ijk-i,fn+i] 

i=0 
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where is a member of the knot vector with n H- A: H- 1 elements 

T = ^ — Ij • • • j ^71— 1j ^n+lj • • • j ^n+A:)j ^0 — ^ ^ 

The basis functions, are given by: 



N^^{u) 



1 for < n < 
0 otherwise 



for A; = 1, and 

N^,k{u) = N,^k-l{n) + 

— 1 ^i+A; 

for A: > 1 and z = 0, 1, . . . ,n. These equations have the following properties [22, p. 168]: 

• ^i,A;(^^) > 0, for ti <U < ti^k 

• = 0, for to <u< ti, and Uj^k <u< tn^k 

• YA=o^ik{y) = 1, for u e [tfc-i, Wi]j and 

• N^k{'^) has continuity at each simple knot ti. 

The first three properties show positivity and partition of unity for the B-Spline basis 
functions. 



p 



r-2 




Figure 3-5: A Cubic B-spline Curve Segment with Control Polygon 

Local control means that a single segment of the B-spline curve is controlled only by the 
nearest n points, and that any control point affects the nearest n spans. In other words, 
changing pi affects the curve in the parameter range U <u < tij^k and the curve at a point 
u where tr < u < is determined completely by the control points • • • ? Pr* 

The convex hull property for B-splines applies locally, so that a span lies within the 
convex hull of the control points that affect it. This provides a tighter convex hull property 
than that of a Bezier curve, as can be seen in Figure 3-6. 
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Figure 3-6: Convex Hull Property for a Cubic B-Spline Curve 



Increasing the multiplicity of a knot reduces the continuity of the curve at that knot. 
Specifically, the curve is {k—p—1) times continuously differentiable at a knot with multiplic- 
ity P k)^ and thus hcis continuity. Therefore, the control polygon will coincide 

with the curve at a knot of multiplicity fc — 1, and a knot with multiplicity k indicates C~^ 
continuity, or a discontinuous curve. 

As you can see in Figure 3-5, the endpoints of the curve may not coincide with the 
control polygon. However, repeating the knots at the end k times will force the endpoints 
to coincide with the control polygon. Thus a knot vector described by 

T • • • ? ^n—l ) ^Ti) • • • > ^n+A:) 

' V ' V ' V ^ 

k equal knots n — A: -f- 1 internal knots k equal knots 

will have to = = . . . = tjt-i 3,nd tn+i = ^n+2 = . . . = tn-fJb ^nd hence the first and last 

control points of the curve coincide with the endpoints of the curve. 

From this discussion, it can be seen that a Bezier curve of order k (degree A: — 1) is a 
B-spline curve with no internal knots and the end knots repeated k times. The knot vector 
is thus 



T (to> tj , . . . , 1 . . . , 

V III ,^..1 .✓ V*" ' '' 

k equal knots k equal knots 
where n-hA:-l-l = 2/;;orn = fc — 1. 

A B-spline surface is a tensor product surface using a B-spline basis. This is represented 
mathematically as 



m n 

R(u,u) = Y,Y^PijNi,k{u)Nj^i{v). 



i=Q j=0 
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Knot Insertion 

A knot can be inserted into a B-spline curve without changing the geometry of the curve. 
The new curve is identical to the old one, with a new basis where 

n n-f-1 

becomes 

t=0 i=0 

over T = [to, ti,..., . . .] over T = [to, ti,..., U, i, tf+i, . . .] 

when a new knot i is inserted between knots ti and The deBoor algorithm is used to 
insert the knot, and is represented as 

Pi = (1 +"iPr^ i>l-k + 2 (3.4) 

where 

ai = ^ 



To evaluate a B-spline curve at a specific parametric value or to split the curve at that 
value, the deBoor algorithm is used to repeatedly insert the same knot until the control 
polygon coincides with the curve at the parametric value uq, or in other words, until the 
multiplicity of the knot at uq = A: — 1. Then, 

P? = Pi 
Pfcl} = r(uo). 

A B-spline curve is continuous in the interior of a span. Inserting a knot does not 
change the curve, so it does not change the continuity. However, if any of the control points 
are moved after knot insertion, the continuity at the knot will become where p is 

the multiplicity of the knot. Figure 3-7 illustrates a single insertion of a knot at parameter 
value uo, resulting in a knot with multiplicity of one. 



1 

0 

UQ-tj 

U+k — 1 



i ^ I — k 1 
Z + 1 ^ 2 
I — + 




Figure 3-7: The deBoor Algorithm 
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The B-spline curve can be subdivided into Bezier segments by knot insertion at each 
internal knot until the multiplicity of each internal knot is the degree k. The above algorithm 
is also known as Boehm^s algorithm [4, 5]. 

Knot Removal 

Knot removal is the reverse process of knot insertion. This thesis uses the knot removal 
algorithm developed by Tiller [42], in which a knot is removed if and only if it leaves the 
curve or surface geometrically and parametrically unchanged. 

To demonstrate the process, this example uses a cubic B-spline curve r(u) given by 
control points (pg, . . . ,Pe) and knot vector (ig? • • * ^lo) where ig = . . . = ^3 = 0, ^4 = ^5 = 
^6 = 1 and ^7 = . . . = ^ig = 2 as shown in Figure 3-8. As the basis functions only guarantee 
continuity at u = 1, the first derivative may or may not be continuous there. Using the 
continuity condition (3.2), the first derivative will be continuous if and only if 



{tj - <4)(P3 - P 2 ) = (<6 - <3)(P4 - Pa)- 



Since t .4 = = u, 



n0_ ^ 

P 3 — 



h — h 



t3_o , tj -u Q 

P4 + — rP2- 



U — h 



Since p^ = p^ and p° = P3, 



P 3 = «3P3 + (1 - «3)P2 «3 = 7 — (3.5) 

17 — 13 

Note that (3.5) is the same as the deBoor algorithm equation (3.4) for inserting a new knot 
at u = 1, although the points and knots are numbered differently in this example. 

A similar reasoning yields that fact that a knot u = I can be removed a second time if 
and only if the second derivative is continuous, yielding 



P2 = «2P2 + (1 - «2)Pi 
P3 = «3P3 + (1 - «3)P2 

^i = -r^—^ * = 2,3 (3.6) 

and a knot it = 1 can be removed a third time if and only if the third derivative is continuous 
such that 



P? =«iPi + (1 -Q:i)Po 
P2 = 02 P 2 + (1 - 02)p? 

P3 = «3P3 + (1 - «3)P2 

= i = 1,2,3. (3.7) 

Note that there are no unknowns in equation (3.5), one unknown, P2, in equation (3.6) and 
two unknowns, pf and P2, in equation (3.7). 

For the knot removal process, first the right hand side of equation (3.5) is computed and 
compared to P3. If they are equal within a given tolerance, the knot and p® are removed. 
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If the first knot removal is successful, equations (3.6) are solved for P 2 and compared: 

- (1 - a2)p\ 

OC2 

P 3 - Q3P3 
1 — as 

If the two values for p| are the same, then the knot and control points P 2 and P 3 are 
removed and control point p| is inserted. 

If the second knot removal is successful, the third step is to solve the first and third 
equations of (3.7) for 




- (1 -ai)pg 

ai 

Pj - Q3P3 

1 - as 

The two values are then substituted into the second equation of (3.7). If the result is within 
tolerance of P 2 , then the knot is removed and control points p^, P 2 and p| are replaced by 
pf and P 2 - 




1 2 2 

P3 P2 P3 




Figure 3-8: Knot Removal 

This can be generalized to apply to any number of removals of any particular knot. For 
the nth removal, there will be a system of n equations with n — 1 unknowns. If n is even, 
two values of the final unknown control point will be calculated and compared. If they 
are within tolerance, the knot removal is successful. If n is odd, all new control points are 
computed and the final two are substituted into the middle equation. If the result is within 
the tolerance, the knot removal is successful. If the nth removal is successful, n control 
points will be replaced by n — 1 control points. 

Knot removal from a surface is performed on the m -f- 1 rows or n -f- 1 columns of control 
points. However, the knot removal is successful only if the knot can be successfully removed 
from each row or column. Therefore, the result must be checked for each row or column 
before any control points are removed. 



34 



NURBS 

The non-uniform rational B-spline (NURBS) curve is represented in a rational form 



where > 0 is a weighting factor and Ni^ki'^) defined over a non-uniform knot vector. 
A NURBS surface patch can be represented as 

where Wij > 0 is a weighting factor. A uniform knot vector is one where the knots are 
evenly spaced. If the knots are not evenly spaced, the knot vector is non-uniform. Thus, 
any knot vector with repeated knots is non-uniform. The NURBS representation of curves 
and surfaces allows the exact representation of figures such as circles, conics, quadratics, 
and surfaces of revolution with rational profiles such as torii. If = 1 for all i or W{j = 1 
for all i,jj this representation is the same as the integral non-uniform B-spline described 
previously. In this thesis, all the B-splines used are integral, but the development could be 
extended to include NURBS. 
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Chapter 4 

Review of Differential Geometry 
Properties of Developable Surfaces 

4.1 Introduction 

This chapter discusses the properties of ruled and developable surfaces. The representations 
and properties discussed herein are taken from classic texts by Struik [41], Kreyszig [26], 
Faux and Pratt [11] and Spivak [40]. 



4.2 Ruled and Developable Surfaces 

Ruled Surfaces 

A ruled surface is a surface generated by the motion of a straight line (a generator or ruling) 
through space [41]. A curve that passes through all the rulings of the surface is called a 
directrix. Any point on the surface can be expressed as 

R(u,i;) = q(u) + v(i{u) (4.1) 

where q(u) is a directrix or base curve of the ruled surface and (^{u) is a unit vector which 
gives the direction of the ruling at each point on the directrix. Alternatively, the surface 
can be represented as a ruling joining corresponding points on two space curves. This is 
represented by 



R(u,u) = (1 — u)r^(u) -f ura(u) 0 < < 1 (4.2) 

where r^(u) and ra(u) are directrices, as shown in Figure 4-1. The two representations are 
identical if 



a(u) = r^(u) and /3(u) = rs(u) - r^(u). (4.3) 

Throughout the thesis it is assumed that the constant u isoparametric lines correspond to 
the generators of the developable surface or, in other words, the straight line rulings are in 
the V direction. 
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directnx, fA(u) 




b2 f 



Figure 4-1: A Ruled Surface 



Developable Surfaces - A Subset of Ruled Surfaces 

Developable surfaces are a subset of ruled surfaces that have a constant tangent plane along 
each generator. Since surface normals are orthogonal to the tangent plane and the tangent 
plane along a generator is constant, all normal vectors along a generator are parallel. This 
is shown in Figure 4-2. 




Figure 4-2: A Developable Surface with Tangent Plane along a Ruling 



Given two points on a single generator that are selected at parameter values (uo,t^i) 
and the surface is developable if the tangent planes at these points coincide. The 

tangent planes are defined by the vectors and R^. Using the representation of a ruled 
surface shown in (4.2), we obtain 



R«(uo,ui) = (1 
'Ru{uq,V2) = (1 - V2) 



drsiuo) 

+ Wl , 

du 


(4.4) 


dreiuo) 

+ V2 , 

du 


(4.5) 
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Kv{uo,vi) = Rv{uq,V 2 ) = VBiuo) - r^(uo). (4.6) 

If the tangent planes at the two points coincide, the vectors R„(uoiVi), ~Ru{uo,V 2 ) and 
rs(uo) ~ ^a{‘^o) must be coplanar. Therefore, the triple product must equal zero: 

|R^,(uo,ui) R„(uo,U 2) (rfl(uo) - rA(uo))l = 0. (4.7) 



Let R„(uo,vi) = R^,l and Ru(uo,U 2 ) = Ru 2 - Then 

R,1 X R,2 = [R1:iR^2 - - [R-SiRu2 - KiK2]S+ " ^111^22]^ (4-8) 



Expanding the first term of (4.8) by employing (4.4) and (4.5) yields 



R y 02 r >2 TyV 

u\^u2 ~ ^ul^u2 



Similarly, 



[(1 - ui)r^ + uir^][(l - V2)t\ 4- U2rB] 

-[(1 - v\)t\ + uir|][(l - V2)t^a + ^ 21 - 0 ] 

(1 - ui)(l - V2)ir^^rA + ~ V2)r%i-B + V2(l - vi)P^r% + viV2r^srB ~ 

[(1 - t;i)(l - V2)r\r\ + vi(l - V2)t\t% + V2{1 - + viV2t\t%] 

(vi — UlU2 —V2 + ViV2)t\t^b + (^2 ~ ^1^2 — Vi + V\V2 )t\t% 



RJiR^2 - K 1 K 2 = (^2 - wi)(r^r^B - 



and 



Therefore, 



~ Ki^l2 = {V2 - vi){r^Ar% - r^r|). 



Rwi X = {V2 - X Tb) 



and thus (4.7) can be reduced to 

(fb(u) - r 4 (u)) X F4(u) • rB(w) = 0 (4.9) 

for any u as long as 7 ^ U 2 . Conversely, any surface that satisfies (4.9) must satisfy (4.7) 
and, therefore, the tangent planes must be the same and the surface must be developable. 
Thus, a surface is developable if and only if (4.9) is satisfied. Substituting (4.3) into (4.9) 
and using the fact that x = 0 yields the equivalent condition 



d X /3 • /3 = 0. 



(4.10) 



Developable Surfaces - Envelopes of Families of Planes 

An alternative representation of a developable surface is as an envelope of a family of planes. 
The concept of an envelope will be described here for a family of curves, then extended to 
planes. First, we represent a curve as r(u,a) where u is the curve parameter and a is the 
family parameter. In other words, each curve in the family is r(u, a^i). There may be a 
curve Fe(a) that is tangent to each curve in the family, as shown in Figure 4-3. This curve 
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would be the envelope of that family of curves [11]. 




Figure 4-3: An Envelope of a Family of Curves 



Similarly, a surface can be the envelope of a family of surfaces. Let a one-parameter 
family of surfaces be represented by G(x,a), where a is the family parameter, and assume 
that the surfaces are twice continuously differentiable. Also assume that consecutive surfaces 
intersect one another. The intersection is 



G(x, a) = 0 G(x, a + /i) = 0 



where |h| is sufficiently small. This may also be represented by 

G(x,g + h)~ G(x,a) _ 



(4.11) 



The limit of (4.11) cis h tends to zero is dG{x,a)/da = 0. The set of points represented by 



G(x, a) = 0 



dG{x, a) 
da 



is the characteristic of the surface. If the characteristics of all the possible surfaces in the 
family form a surface, that surface is called the envelope. If the envelope exists, then at 
every point of the characteristic of a surface in the family, the tangent planes of the surface 
and of the envelope coincide. 

A characteristic point is determined by the intersection of the characteristic with another 
surface in the family. These points of intersection are determined by 



G(x,a)=0 ^<j(x, Q^) _ Q G{x,a + k)=0. 

oa 

where |A:| is sufficiently small. Using the Taylor series expansion 

t\ ru \ ,dG{x,a) d‘^G{x,a 9k) 

G(x, g + A;) — G(x, a) -i- k h s~9 1 0 <C 9 <C 1 , 



da 



da^ 
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the third equation can be rewritten as 



^G{y.,ot + 6k) 
da^ 



Then, taking the limit as k approaches zero, the characteristic points of the surface are 



G(x, O') = 0 



dGj-x., g) ^ a^G(x,a) 

da da^ 



If the characteristic points exist and if they form a curve as a varies, the curve lies on the 
envelope of the family and is called the edge of regression of the envelope. If the envelope 
and the edge of regression exist, then, at every point of the edge of regression, the tangents 
of the edge of regression and the corresponding characteristics coincide. The characteristics 
then form two sheets corresponding to the tangents of the edge of regression in the positive 
and negative directions. These sheets form a cusp at the edge of regression, as shown in 
Figure 4-4. 




Figure 4-4: Edge of Regression 



by 



Let us now consider a family of planes and their derivatives which can be represented 



G{x, a) 

da 

^G 

da^ 



a(o:) • X + c(a) = 0 
dsi dc ^ 

d^a. d^c 



where a(a) is the unit normal vector to the corresponding plane. When all planes are paral- 
lel, no intersection exists and no characteristic curves are formed. When all planes intersect 
in a single line forming a pencil, the envelope is a line instead of a surface. Therefore, these 
two cases are excluded [26]. When the envelope exists, it is formed of the characteristics of 
the family of planes. Since planes intersect in straight lines, the envelope is a ruled surface. 
Furthermore, the tangent plane along a characteristic is constant since this is an envelope 
of planes, and hence the surface is developable. 
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When the vectors a(a) are all parallel to a single plane P, then all the planes in the 
family are parallel to a vector normal to the plane P. Thus, the envelope of the family of 
planes forms a cylinder (except when they form a pencil which has been excluded). When 
all the characteristic points coincide, all the characteristics intersect in a single point and 
the surface is a cone. When the characteristic points form a curve, the curve is the edge 
of regression and the surface is the tangent surface to the curve. Therefore, all developable 
surfaces can be represented as planes, cylinders, cones and tangent surfaces, or combinations 



thereof [41]. 

This description of developable surfaces 
curves and developable surfaces [41, p.72j: 

Curves 

2 points determine a line. 

3 points determine a plane. 

2 consecutive points on a 

curve determine a tangent line. 

The curve is the envelope 
of tangents. 

Gaussian Curvature 
Using equation (4.1) 



leads to the following duality between space 

Developable Surfaces 

2 planes determine a line. 

3 planes determine a point. 

2 consecutive planes of the 

family of planes determine 
a characteristic line. 

The developable surface is 

generated by characteristic lines. 



R(w, v) = a(u) + v0{u)^ 



the first and second derivatives of a surface can be represented as 



Ru(u,u) = 
R„(u,u) = 
Ruu (^) y) ~ 

Also recall from (2.9) that 



d(u) -h v/3{u) 


(4.12) 


(3{u) 


(4.13) 


a{u) + v(3{u) 


(4.14) 


0{u) and 


(4.15) 


0. 


(4.16) 


X > 0 





EG-F^ = |R„ 
for a regular surface. The surface normal is given by (2.7) 



_ X R„ 

jRu ^ Ri;| 

_ (d(u) + v0{u)) X (3{u) 

VEG- F2 

_ 6l{u) X /3(ii) + v{0{u) X /3(u)) 

VEG -F^ 

d(u) X (3{u) 

VEG - F2 ■ 
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Recall from (2.23) that Gaussian curvature is represented as 



j. _LN - A/2 

^ “ EG -F^' 

where, from equations (2.13) and (2.12) 

N = Rt,„ • N and 
M = Rt,„-N. 



Thus, from (4.16), N = 0 and the Gaussian curvature can be rewritten as 



K 



m2 

EG-F^ 

(Ru. • N)2 
EG-F‘^ 

i${u) • a{u) X /3(u))2 
{EG - F2)2 
(gx/3-/3)2 
{EG - F2)2 • 



(4.18) 



Thus, from condition (4.10), zero Gaussian curvature is a sufficient and necessary condition 
for a surface to be developable. 

Since the Gaussian curvature of a developable surface is zero everywhere, the maximum 
and minimum principal curvatures (2.21) and (2.22) of a developable surface can be written 
as 



^max — H \H\^ f^min — H 
These principal curvatures reduce to 



max — 2iif, 


^min ~ 0 


when 


H>0, 


(4.19) 


^max — Oj 


^min — 0 


when 


H = 0 and 


(4.20) 


^max “ Oj 


^min — 


when 


o 

V 


(4.21) 



It is clear from equations (4.19) to (4.21) that at least one of the principal curvatures is 
always zero on a developable surface, which agrees with the fact that the Gaussian curvature 
is zero everywhere (see equation (4.18)). «max and ^min from (4.19) and (4.21) respectively 
are each termed the nonzero principal curvature^ k*, where 

k" = 2H, (4.22) 



Mapping 

Isometric mapping is defined as a mapping that preserves arc lengths between the mapped 
surfaces. Two surfaces that can be isometrically mapped to one another are called isomet- 
rics. Since arc lengths are preserved, the first fundamental forms of isometrics are the same. 
Since the Gaussian curvature of a surface and the geodesic curvature of a curve on a sur- 
face depend only upon the coefficients of the first fundamental form, corresponding points 
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on two isometric surfaces have the same Gaussian curvature and corresponding curves on 
isometric surfaces have the same geodesic curvature at corresponding points [26]. 

A conformal mapping is one in which the angle between every pair of intersecting arcs 
is the same on both surfaces. This occurs if the coefficients of the first fundamental form on 
the original surface are proportional to those on the inverse surface [26]. Therefore, every 
isometric surface is conformal. 

Since a plane has zero Gaussian curvature, only a surface with zero Gaussian curvature 
can be mapped to it. Therefore, only a developable surface can be isometrically mapped to 
a plane [26, p.l89]. 

Since the tangent planes along any ruling are constant, developable surfaces can be un- 
rolled, or developed, isometrically into a plane without stretching or tearing by successively 
placing each ruling on the plane. One can think of a developable surface as a paper with 
folds. As the number of folds approaches infinity, the surface becomes smooth. The surface 
can then be unfolded again into a plane. This is shown in Figure 4-5. 




Figure 4-5: A Surface with Folds and a Developable Surface 

A geodesic represents the curve with the shortest length between two points on a sur- 
face. Therefore, geodesics on a plane are straight lines. Since developable surfaces map 
isometrically to planes, geodesics on developable surfaces map to straight lines on a plane 
and straight lines on a plane map to geodesics on a developable surface. 

Developable Surface Properties 

The following statements are equivalent necessary and sufficient conditions for a surface to 
be developable. Each has been described above. Thorough proofs can be found in Struik [41] 
and Kreyszig [26]. 

1. Developable surfaces possess the same tangent plane at all points of the same gener- 
ator. 

2. The normal vectors on a developable surface along a ruling are parallel. 

3. A developable surface is the envelope of a one-parameter family of planes. 

4. Gaussian curvature is zero. 
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5. The mapping of a developable surface onto a plane is isometric. 

6. Geodesics on a developable surface map to straight lines on a plane, and straight lines 
on a plane map into geodesics on a developable surface. 
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Chapter 5 



Design With B-Spline Developable 
Surfaces 

5.1 Introduction and Literature Review 

In the current literature, there are three main approaches to representing developable sur- 
faces in Bezier form. These approaches are described below. This chapter explores user- 
friendly methods for the design of developable surfaces. Section 5.2 describes the design of 
developable B-spline strip surfaces that lie between parallel planes. Section 5.3 describes 
the design of developable B-spline surfaces with general three-dimensional directrices. The 
development of a developable surface onto a plane is described in Section 5.4. This is 
an efficient implementation of the method in the current literature, also described below. 
Section 5.5 provides some engineering examples. 



5.1.1 Developable Bezier Strips 

Aumann [1] provides the basis for a user friendly and computationally inexpensive method 
of designing developable surfaces. He represents developable surfaces in terms of two Bezier 
curves (directrices) and rulings between pairs of points from each curve, as discussed in 
Section 4.2, but restricts the two directrices (r/i(n), rB(u)) to lie on parallel planes so that 
the tangent vectors r>i(u) and vb{'^) are parallel, i.e. where p{u) is a 

linear function of u. This condition automatically satisfies the condition (4.9). 

His design philosophy is to consider r^(u) as a design curve, which means its degree, knot 
vector and control points are specified. In addition, the two end points r^(0) and tb{1) of 
the resulting directrix are provided as boundary conditions. The knots and internal control 
points of the second directrix vb{u) are then computed using a simple system of linear 
equations. He investigates the case when the design cmve is a cubic Bezier curve and p{u) 
is a linear function of u. Hence, the resulting directrix is quartic. Figure 5-1 illustrates the 
terminology. 

Frey and Bindschadler [14] extend the work of Aumann [1] by generalizing the degree of 
the directrices and by considering the case where p{u) is quadratic. Aumann [2] further ex- 
tends his work so v/i{u) and vb{u) have the same degree. All three papers discuss regularity 
conditions of the developable surfaces. 
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design curve, r^Cu) 




b2* 



Figure 5-1: A Developable Surface 

5.1.2 Duality Between Points and Planes 

The second approach, introduced by Bodduluri and Ravani [3], is based on the represen- 
tation of developable surfaces as an envelope of a one-parameter set of planes as described 
in Section 4.2. They express the surfaces in terms of plane geometry using the concept of 
duality between points and planes in 3D projective space. They interpret the developable 
surface as the set of its tangent planes represented in the dual form. Their method is math- 
ematically elegant; however, it is somewhat less user friendly and requires a conversion to 
a standard tensor product form [37]. 

Pottmann and Farin [37] generalize the results of Bodduluri and Ravani [3] by providing 
algorithms for converting the dual representation of developables to standard tensor product 
form. Hoschek and Pottmann [23] derive algorithms for interpolation or approximation 
with developable B-spline surfaces. Given a developable surface patch D, they pick out 
an appropriate number of tangent planes Yj of a developable surface patch D, and either 
interpolate or approximate this data with a developable B-spline surface. 

5.1.3 Nonlinear Representation 

The third approach is due to Lang and Roschel [29]. They obtain conditions for the control 
nets and weights of rational Bezier surfaces to be developable, which leads to a complicated 
nonlinear system of conditions. Although their work provides a useful criterion for checking 
the developability of a given ruled surface patch, it is difficult to design a developable surface 
with this method. 

5.1.4 Development 

Faux and Pratt [11] provided a method for the development of a curve on a developable 
surface onto a plane based on the fact that a developable surface maps isometrically onto a 
plane. This method is also described by Gurunathan and Dhande [16]. Another method is 
described by Clements and Leon [7], in which a geodesic on the surface is mapped to the 
plane and points on the plane are determined in reference to this geodesic. This method. 
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however, requires the mapping of the geodesic and then the calculation of the location of 
each edge, and is therefore slower than the other method [11]. 



5.2 B-Spline Developable Surfaces with Directrices in Paral- 
lel Planes 

This section further extends Aumann’s work [1] such that n Bezier developable surfaces, 
whose two directrices lie in parallel planes, are strung together by joining them along their 
end rulings. In this manner, they can be represented as a single B-spline surface with 
continuity. 

From Section 4.2, a surface is developable if and only if the tangent plane is constant 
along a generator. If the two directrices r>i(u) and r^(u) are restricted to lie in parallel 
planes, then the tangent vectors to r^(u) and rB(^) must be parallel at every u for the 
surface to be developable. Therefore, condition (4.9) is equivalent to 

FB ( ti ) = p{u)va{u) (5.1) 

where p{u) is a scalar function of u. In this thesis, p{u) is selected to be linear, expressed 
as 



p{u) = po(l -u) + piu. 

The design philosophy is to consider r^(u) as a design curve [1, 14], which means its 
degree, knot vector and control points are specified. Also, the two end points r 0 (O) and 
fb( 1) are provided as boundary conditions. The knots and internal control points of the 
second directrix, fb(u), are then computed. Since multiplication of B-spline functions is 
complex and computationally expensive, the Bezier segments are first extracted from the 
B-spline curve by knot insertion using Boehm’s algorithm [ 22 ] described in Section 3.3. 

Figure 5-2(a) shows the design curve r^(u) which is a planar cubic B-spline curve with 
knot vector T=(0 0005 |llll) and control points (0, 2, 0), (0.8, 0.5, 0), (1.6, 1.5, 
0), (2.4, 1.5, 0), (2.6, -0.5, 0) and (4, 0, 0). Figure 5-2(b) shows three Bezier segments 
after knot insertion with the two endpoints of re(u), which are supplied by the user to be 
tb( 0) = (—1, 3, —3) and r^(l) = (2, 1.3, —3). Note that Fyi(u) is in the z = 0 plane and the 
F^ endpoints are in the z = 3 plane, parallel to the plane of r^. 

Let F^(u) be a degree n curve with m Bezier segments. Since p{u) is linear, the degree 
of fb(u) is n + 1. We denote the jth Bezier segments as 

n 

j = 0,...,m- 1 and 

t=0 

n+1 

r^(u) = J = 0, . . . , 771 — 1. (5*2) 

i-0 

The control points a and b are two-dimensional vectors since they have been restricted to 
lie in parallel z planes. Substituting (5.2) into equation (5.1) yields 

n 

(n+l)EO,( 

i=0 
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(a) Design Curve r^(u) with Controi Potygon and Knots (b) Design Curve after Knot Insertion with Bezier Segments 

and Endpoints of Tg(u) Shown 




Figure 5-2: Design of Developable B-Spline Strip Surface 



i=0 

n — \ 

= Tl ^ ^ (PO ( 1 u) -f- (5.3) 



1=0 



for j = 0, . . . ,m — 1. Using the definition of a Bernstein polynomial (3.1), 



(1 -'u)Si,„_i(u) = (1 -u)t 



in 1)1 (1 _ 



i\{n — 1 — i)! 
n — 1 (n)! 



n 2 !(n — i)! 



(1 - 



n 



(5.4) 
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and 



7-, / \ (^ ~ !)• /I \ 

i + 1 n! 



n-l-i^i 



n (i + 1)!(^ ~ (* + !)• 
i + 1 



j(l 



n 



■Bt+l,n{u). 



(5.5) 



Substituting (5.4) and (5.5) into (5.3) yields 



n 



(n + 1) y^(b^-(n+i)+i+i 
i=o 



t^j(n+l)+i)'®i,n(^) 
n— 1 

n ^ [po(l - u)Bi^n-i{u) + piuBi^n-i{u)] (a^n+i+i - aj„+i) 

1=0 

n-1 

^ [po(n - i)Bx,n(u) + Pi{i + l)Si+i,n(u)] (aj^+i+i - ajn+i) 
i=0 

Po^(ajn+i aj>i)5o,ri(^) 

n — 1 

+ [po(^ “ 0(^jn+i+i ” aj>i_|_2) + pii(ajyi-)_i — a^Ti-i-i—i)] Bi^ji{u) 

f=l 

4"Pl^(ajn+n ^jri'^n — l)^nyn{'^)^ 



Since Bernstein bctsis polynomials are linearly independent, each set of coefficients must be 
equal, yielding the following system of n — 1 equations [14] for the jth Bezier segment 

^Poi^jn-\-i ajn) (n + 1) (5.6) 

(n — i)pQ{a.jn-\-i-\.i — Q,jn-\-i) + ip{{Q^jn-\-i ~~ aj^+i-l) = (n + 1) — bj(„_|_i)_|_j) 

^Pl (^(7 + l)n ^(j+l)n— l) (t2 + 1) ), 

where z = 1, . . . , n — 1. 

Now consider joining two consecutive Bezier patches with continuity. First, 
continuity requires that the first derivatives along the connecting edges of both patches 
must be equal. Since the connecting edge is a ruling and is therefore linear, this requirement 
reduces to 



^a(I) = (5-7) 

4(1) = 4"'(0)- (5-8) 

From equation (5.1) we have 

4(1) = />^'(i)4(i) 

4‘"‘(0) = p^+i(0)4"'(0)- (5-9) 

Substituting equations (5.9) into equation (5.8) and using equation (5.7) yields 

(5.10) 
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continuity requires that the second derivatives along the connecting edges must be 
equal. Again, the connecting edge is linear, which reduces this requirement to 

i^^(l) = i^+‘(0) (5.11) 

(5.12) 

Taking the derivative of equation (5.1) yields 

tb[u) = p{u)va{u) + p{u)TA{y) (5.13) 

where 

p{u) = Pi- Po. 

Evaluating equation (5.13) at the connecting edge yields 

i^B'"‘(0) = + (5.14) 

Substituting equations (5.14) into equation (5.12) and using equation (5.11) yields 

2p"o^' = + (5.15) 



Each Bezier patch has n -f 1 control points along r^(u) and n + 2 control points along 
tb(u), where n is the degree of the original B-spline design curve. The system of equations 
for a single patch derived in equation (5.6) can be solved for the coefficients of the second 
directrix 






^ ^Poi^jn+l 3-jn) 



t>j(n+l)4-i+l 



bO + l)(n+l) 



bj(n+l)+i ^ ^jn+i) 

i 

i = 1, . . . , n — 1 

Tl 

^(i+l)(n+l)-l + „ + “ ®0 + l)n-l) (5.16) 



where j = 0, . . . , m — 1 for a system of m patches. Adding equations (5.16) for m patches 
yields a vector equation that can be solved for terms of the first and last 

points of the second directrix: 



where 



bm(n+i)-bo = PoVo+pV" ^Vi 



1 \ 

Vo = 2^ a, + ^ ain - 



ao 



and 



m—l 



Vl = 



^mn 



m(n + 1) 



E^+ E 



i=l 



(5.17) 
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Combining equations (5.10) and (5.15), we can express and in terms of the two 
end coefficients pg follows: 

/^i~^ j = 0,...,m-l. (5.18) 

Therefore, given a design directrix, F 4 (u), and the endpoints of the second directrix, r^(0) 
and rB(l), all the interior control points of the second directrix can be obtained succes- 
sively using equations (5.16). After the second directrix is calculated, the first directrix is 
degree elevated so the two curves may be combined into a B-spline strip. Any knots that 
may be removed from the surface while leaving the curve geometrically and parametrically 
unchanged are removed using the knot removal technique [42] described in Section 3.2. 

Figure 5-2(c) shows the computed rB(u) together with r^(u) in Bezier form and the two 
end rulings. Figure 5-2(d) shows the final B-spline developable surface after knot removal. 

Developable surfaces become irregular when |R„ x R^| vanishes, which leads to the 
regularity condition for a strip surface [14] 

p{u) >0, 0 < u < 1. 



From a user’s point of view, once he/she determines one of the endpoints, it is convenient 
to know the range that the other end point can take. From (5.18) it can be seen that 
if both Po and are positive, then all p coefficients are positive and the surface is 

regular. The coefficients p can be determined by solving the vector equation (5.17). The 
resulting equations for pg and p^~^ must be positive for the surface to be regular, yielding 
the following restrictions: 



v^v\-vfvS 






(5.19) 



and 



_1 _ (^m(n+l) - ^o)^0 - (^m(n+l) ~ ^o)^0 



Pi = 






v;^v{ - vrv, 






>0 



1 ^0 



(5.20) 



where the superscripts denote the x and y components of the vectors. Given one endpoint, 
a regular surface may be formed by selecting the second endpoint to lie within the range of 
values determined by equations (5.19) and (5.20). 

It is also interesting to observe [14] that when the two vectors joining two end points of 
each directrix are parallel and their magnitudes are multiples of a constant, the resulting 
developable surface is a generalized cone. When the constant is equal to one, the surface is 
a cylindrical surface. In other words, if p{u) is a constant not equal to one, the surface is a 
generalized cone, and if p{u) = 1, the surface is a cylinder. 



Algorithm 

The following steps describe the algorithm for the solution method. 

1. The user provides the design curve r^(u) and the endpoint re(0). The program 
provides a range of values for Fb( 1) that ensure a regular surface, and the user supplies 
a specific value for rB(l). 
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2. Bezier segments are extracted from r^('u) by knot insertion. 



3. The system of equations (5.16) is solved for the internal control points of tb{u), 

4. Once the resultant directrix ra(u) is computed, the developable Bezier surface is 
obtained by degree raising r>i(u) and substituting in equation (4.2). 

5. Knots that can be removed while leaving the surface unchanged are removed using 
the method described in Section 3.2. 



5.3 Developable Surfaces with 3D Directrices 

5.3.1 Developable (n, 1) Bezier Surfaces 

In Section 5.2, we restricted the two directrices to lie on parallel planes. In such case, r>i(?i) 
and re(u) must be parallel to satisfy condition (4.9). In this section, the directrices are not 
constrained to lie on parallel planes, so equation (5.1), = p(u)r^, no longer applies. 

The design philosophy here is the same as that of the previous section. Namely, r^(u) 
is designated the design curve with its degree, knot vector and control points specified. The 
two end points r^(0) and rfi(l) are provided as boundary conditions and the knots and 
internal control points of the second directrix, tb{u) are computed. 

We begin this development by assuming the design curve r^(u) is a Bezier curve of 
degree n, 



^/l(^) — ^ ^ (^) ? 0 ^ ^ Ij 

i=0 

where = (af,af,a^)^ are control points. The resultant curve rB(u) is also a Bezier curve 
of degree n. In Section 5.3.2 we investigate the case when the curves r^(ii) and r^(ii) are 
B-spline curves. If hi = (5f,6f,6f)^, i = 0, ...,n, are the control points of the degree n 
Bezier curve vb{u)^ and ai, bj and Cj are 



pi 


II 

II 

-f 

1 


i = 0,.. 


. , n — 1 


bi 


= (6f,6V,6n"’ = bi+i-bi 


II 

o 


. , n — 1 


Ci 


= (cf,cf,c?)^ = bi -aj 


1 = 0,.. 


.,n. 



then equation (4.9) can be expressed as 

CiSi^^('u)'^ X ^Ti ^ ^ X (u)^ • (u)^ =0. (5.21) 



Ki=0 



1=0 



Since equation (5.21) is a scalar polynomial equation of degree 3n — 2, we can rewrite it 
(upon division by n^) in the form 



3n-2 

^kBk,3n-2{'^) = 0 

k=0 



(5.22) 



54 



where 



and 



min(2n— 1,/c) /2n— Iwn— 1\ 

Bk= E + 

i=max(0,fc— n+l) 



(^-2) 



/n\ /n-l\ 

^i-j * 
(2n-l) 

mv-]) 

/2n-lx 

j=max(0,i— n+l) ^ i ' 

min{n,i) fn\ /n— 1\ 

j=max(0,i— n-l-1) v i / 



min(n,t) /n\ 

Pf= E 

j=max(0,i— n-f*l) 

min(n,i) /'n — 1\ 



(5.23) 



using the algorithms for multiplication of polynomials in Bernstein form described by 
Farouki and Rajan [10]. Since equation (5.22) is valid for all u and the Bernstein basis 
functions are linearly independent, must be zero for all A:, 0 < A; < 3n — 2. This leads us 
to a system of 3n — 1 polynomial equations with 3(n— 1) unknowns, i.e. (i = 1, . . . , n — 1) 
where the two boundary points, bo and b„, are given. Therefore, there are two more equa- 
tions than unknowns and the system is overdetermined. The system of equations was solved 
using a Numerical Algorithms Group program [35] for minimization. 



Solution Method Using Unconstrained Optimization 

One way to solve the overdetermined system is to seek a least squares solution to E(x) = 0 
where 






and 



E(x) = (Eq{x),Ei{x),... ,E3n-2{x)f- 

This can be obtained by minimizing the function 

-j -j 1 3 ti — 2 

F(x) = - (e2(x) + £;2(x) + . . . + El_,{x)) = 2 E^(x)E(x) = i ^ EU^)- (5.24) 

1=0 

The problem can then be reduced to searching for zeros of the gradient vector field VF(x), 
i.e. 



VF(x) = 0. (5.25) 

The Newton- Raphson method can be used to iteratively solve the nonlinear system (5.25). 
The Taylor series expansion about the current point is 

VF(x -I- 5x) w g(x) -I- [H(x)](5x (5.26) 
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where g(x) and [H(x)] are the gradient vector and the Hessian matrix of the objective 
function F(x). Using (5.26), the step direction Jx is found from 

[H(x)]^x = -g(x). (5.27) 

The value for the subsequent iteration, x is 

X = X + 5x. 

Using the special properties of a least squares function such as equation (5.24), the deriva- 
tives of F(x) are represented as 

g(x) = VF(x) = [A(x)]E(x) 



and 



3n-2 

[H(x)l = v2F(x) = (A(x)][A(x)]^ + Yl Ei{K)V^Ei{K) 

z=0 

where 

[A(x)] = [VFo,VFi,...,VF3„_2] 

is the (3n — 3) x (3n — 1) Jacobian matrix, the columns of which are the first derivative 

dE 

vectors VEi of the components of E {[Aij] = *^). 

In any minimization process, the initial approximation for x greatly affects the resulting 
solution. Here we use a linear interpolation to obtain the initial values for b|, where 

Z TX — Z 

bj = aj + -(bo - ao) H (b„ - a„), z = 1, . . . , n - 1. 

n n 

Simple Bounds 

In many applications, it is helpful to introduce simple bounds on the variables [15] such 
that li < Xi < where Ij and are lower and upper bounds, respectively. This is a 
constrained optimization problem; however, it has a particularly simple form. A solution 
X* is termed feasible if it satisfies all constraints. A step Jx is in a feasible direction if it 
does not cause x to exceed any of the bounds. 

The following discussion treats a generic iteration, so the iteration number has been left 
out for simplicity. At any iteration, some of the variables will be equal to their upper or 
lower bounds. These variables Xpx, are fixed on their bounds and are consequently termed 
fixed variables. The associated bounds are placed in the active set The remaining free 
variables XpR are allowed to vary within the bounds. The search direction 5x is chosen to 
change only the free variables. Therefore, Jxpx = 0 and equation (5.27) reduces to 



[Hfr(x)]^Xfr = -gFR(x). 

With <5 xfr we update x for the subsequent iteration as 



(5.28) 



If (JxpR forces XpR to violate one of the bounds, that bound is added to the active set for 
the subsequent iteration. 

At each iteration, the active set must be analyzed to determine if the constraints in the 
active set should remain active. A solution x* is a local minimum only if there is no feasible 
step Sx that results in a lower value. Thus, if 

(Jx)^g(x) (5.29) 

is non-zero for any 5x, then x* is not a local minimum. Then (5.29) can be represented as 

g(x).= [A]^A 



where [A] is the set of constraint coefficients and A are Lagrange multipliers. For simple 
bounds, [A] is composed of 1 for /, and —1 for Uj, so the Lagrange multipliers can be 
estimated as 



f gi(x) for Xi = li 
y -gi(x) for Xi = Ui 



(5.30) 



Thus, if a Lagrange multiplier for an active simple bound is negative, the corresponding 
bound should be removed from the active set for the subsequent iteration. 




Figure 5-3: Effect of Simple Bounds on Developable Surfaces 



Many solutions are possible for a single set of equations. For example. Figure 5-3 shows 
two equally viable solutions to the same set of equations. Both examples have design curves 
with control points (0, 0, 0), (1.8, 3, -0.7), (3.3, -2, 1.5) and (4, 0.5, 0) and endpoints (0, 0, 
3) and (4, 0.5, 3) for the second directrix. Figure 5-3(a) is unconstrained, yielding interior 
control points for vb of (1.8, 3, 3) and (3.3, -2, 3). Figure 5-3(b) is unconstrained in the x 
and y directions, but has a lower bound equal to 5 set in the z direction. This produces a 
resultant curve with interior control points (1.8, 3, 5) and (3.3, -2, 5). Since both cylinders 
are developable surfaces with continuity that include the endpoints specified by the user, 
they are equally valid solutions to the equations. Both surfaces are exactly developable with 
zero Gaussian curvature. Although this example is rather trivial, it demonstrates the control 
available to the user to affect the shape of the resulting surface by specifying bounds on the 
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solution. Another example is given in Section 5.3.2. 



Generalized Elimination Method 

From equation (5.23) it can be seen that Eq{x) = 0 reduces to a linear equation, 

(bo - ao) X (ai - ao) • (bi - bo) = 0, (5.31) 

where bi is the only unknown. This condition forces the first four control points, namely 
ao, ai, bo and bi, to lie on a plane, and has the following geometric interpretation. The 
two Bezier curves r^(u) and r^(u) pass through the points ao and bo and their tangents 
are in the directions of the vectors aoai and bobi. Since a developable surface hcis the same 
tangent plane at all points along a ruling and these tangents are on the same ruling, u = 0, 
they must be coplanar. Therefore, ao, ai, bo and bi must lie on the same plane. 

Similarly, the last equation, £' 3 „_ 2 (x) = 0, is linear in b„_i: 

(bn — an) X (an an— i) • (bn ~ bn— i) = 0. (5.32) 

The linear equations (5.31) and (5.32) can be used to reduce the dimensionality of the 
optimization by two. The new objective function can be introduced without terms Eq{x) 
and £’3n-2(^). i-e- 



F(x) = F^(x) + FICx) + • • • + F 3 „_ 3 (x). 

The following discusses the generalized elimination method [13, 15], in which the number of 
variables in the optimization is reduced by the number of independent linear constraints. 
In general, if we denote the number of variables in the optimization as N (i.e. x E R^) 
and the number of linear equations as M, the set of linear equations can be written as 

[Ajx = b 

in which [A] is M x N and b E In other words, the ith row of [A] contains the 

coefficients corresponding to the zth linear equation, and is therefore the gradient of that 
equation. 

This reduction in dimensionality can be described formally in terms of two subspaces, 
which can be described as 

• the M-dimensional subspace defined by the rows of [A]. 

• the N — M complementary subspace of vectors orthogonal to the rows of [A]. 

Let [Y] denote any N x M matrix whose columns form a basis for the range space of [A]^. 
In fact, [Y] may be taken cts [A]^. Let [Z] denote a,n N x {N — M) matrix whose columns 
form a bctsis for the set of vectors orthogonal to the rows of [A], so that [A][Z] = [0]. 
The matrix [Z] has linearly independent columns Zi, Z 2 ,..., zat_m which are in the null 
space and act as bctsis vectors for the null space. Since [Y] and [Z] define complementary 
subspaces, every AT-vector x has a unique expansion as a linear combination of the columns 
of [Y] and [Z] 

X = [Y]xy -h [Z]x^, (5.33) 
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where M- vector Xy is the range-space portion of x and {N — M)-vector x^ is the null-space 
portion of x. By using the chain rule on the derivative of equation (5.33), it follows that 

Vx. = [Z]^Vx. 



Thus, 

V^,F(x + [Z]x) = [Z]^V^F(x) = [Z]^g(x) 

is the reduced gradient vector and 

VlF{x + [Z]x) = [Z]^v2f(x)[Z] = [Z]^[H(x)][Z] 

is the reduced Hessian matrix. Given a fecisible initial point xq we solve a reduced size 
{N — M) X {N — M) linear system 

|Z|’'[H{x)|lZ|x, = -lZfg(x) 

for x^, where the subsequent iteration is 

X = X + [Z]x2. 



5.3.2 B-Spline Developable Surfaces 

In this section we merge m Bezier patches together by joining them along their end rulings 
with continuity such that the resulting developable surface can be represented by a single 
B-spline surface. The design procediure is the same as that of the Bezier case, except that 
the design curve is a B-spline curve. Since multiplication of B-spline functions is complex 
and computationally expensive, we first extract the Bezier segments from the B-spline curve 
by knot insertion [5] as described in Section 3.3. We then evaluate equation (5.23) for each 
Bezier segment and add the resulting equations to the objective function after squaring 
them. If equations Eq and £^m(3n-2) ^^e used in the general elimination method, they are 
not included in the objective function. 

Recall from Section 3.2 that, for a degree (n, 1) B-spline surface with m Bezier segments, 
the and continuity conditions at the joints are given by 

^i(bfii-|-i byii) = — byii_i), 1 ^ z ^ m (5.34) 



and 



bjli— 1 , (byii— 1 \yni^2) byij-l-1 + , (byii-j-l byii-|-2)* 

fli 



(5.35) 



Equations (5.34) and (5.35) reduce to six scalar linear equations for each junction. Using 
these linear equations in the generalized elimination method described in Section 5.3.1 
ensures that the surfaces will be exactly continuous; otherwise they may have some 
residual value after optimization that would cause discontinuity. For a degree (n, 1) B-spline 
patch with m Bezier segments, the number of unknowns and linear equations is 3(mn — 1) 
and 6(m -1)4-2 respectively, so the reduced number of unknowns is 3m(n — 2) 4- 1. 

For a B-spline surface, the initial approximation for x is calculated in the same manner 
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as for the Bezier case, expressed as 

2 TTXTl — 2 

= H (bo-ao)H (b,„„ - z = 1, . . . ,mn - 1. (5.36) 

mn mn 

Once the control points of the second directrix are computed, any knots that may be 
removed from the curve while leaving the curve geometrically and parametrically unchanged 
are removed using the knot removal technique described in Section 3.3. 

Figure 5-4(a) shows the design curve r^(u) which is a cubic B-spline curve with control 
■points (0, 4, 0), (1, 2, 1), (2, 4, 1), (3, 1, 0) and (4, 2, 0) and knot vector (0, 0, 0, 0, .5, 1, 
1, 1, 1). Figure 5-4(b) shows the design curve after knot insertion with the resulting two 
Bezier segments and the selected endpoints r^(0) and r 5 (l), which are (0, 4, 3) and (8, 0, 3). 
When this problem is treated as an unconstrained optimization, the answer converges to the 
trivial solution where the internal control points of the new curve are equal to the internal 
control points of the design curve as shown in Figure 5-4(c). To avoid this situation, lower 
bounds are set on the z variables where li = 3.0, 2 = 2, 5, 8, 11 and 14, and the resulting 
internal control points are (2, 0, 5), (4, 4, 5) and (6, -2, 3). The resulting surface is shown 
in Figure 5-4(d) with two Bezier patches and the control mesh delineated. Figure 5-4(e) 
shows the final developable B-spline surface after knot removal. 

Aumann [2] proves that a degree (3-1) Bezier surface whose directrices lie in parallel 
planes is developable if and only if one of the following conditions holds: 

1. The rulings R(0, u) and R(l,'y) are intersecting or parallel. 

2. The coefficient of the cubic term of the design curve is zero. 

3. The design curve contains exactly one singular point. 

In constructing examples for this thesis, it was found that skew end rulings induce non-zero 
Gaussian curvature into (3-1) Bezier surfaces whose directrices are 3D space curves cts well. 

Algorithm 

The following steps describe the algorithm for the solution method. 

1. The design curve, r^(tx) and the endpoints r5(0) and F 5 ( 1 ) are provided by the user 
along with any constraints such as simple bounds or constraining planes on the interior 
control points of r 5 (u). 

2. Bezier segments are subdivided out of the B-spline curve by knot insertion as described 
in Section 3.3. 

3. Initial values for x = (6q ? • • • » ^mn) calculated using a linear inter- 

polation method (5.36). 

4. The objective function is optimized in two steps. First, beginning with the initial 
values, an initial feaisible point is calculated that meets all constraints and bounds. 
Then, an optimal solution is found following the method described in Section 5.3.1, 
using the following steps: 

• If X satisfies the objective function F{x) = 0 and the Lagrangian multiplier 
estimates are non-negative for all active bounds, the iteration process is complete. 
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(a) Design Curve r a(u) with Control Polygon and Knots 



(b) Design Curve after Knot Insertion with Bezier Segments 
and Endpoints of r b(u) Shown 




(c) Resulting Surface Without Simple Bounds 




Figure 5-4: Design of Developable B-Spline Surfaces 
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(Note that the objective function is a sum of squares, so the minimum possible 
sum is zero.) 

• A feasible step direction Sx is found using (5.28). 

• If Sx causes x to exceed any of the bounds, those bounds are added to the active 
set. 

• The Lagrangian multiplier estimates are calculated for all bounds in the active 
set using (5.30). If any are negative, the associated bound is removed from the 
active set. 

• X is set to x + Sx and the next iteration is begun. 

5. Once the resultant directrix r^(u) is computed, the developable B-spline surface is 
found using equation (4.2). 

6. Knots that can be removed while leaving the surface unchanged are removed using 
the method described in Section 3.2. 

7. The Gaussian curvature bound is calculated using the method described in Sec- 
tion 5.3.3 to gauge the quality of the resulting surface. 

8. The surface is developed onto a plane using the technique described in Section 5.4. 
5.3.3 Accuracy of Resulting Developable Surfaces 

Since the optimization results do not always precisely satisfy F(x) = 0, we need a measure 
to determine how close the resulting surface is to an exactly developable surface. Recall 
from Section 4.2 that a necessary and sufficient condition for a surface to be developable is 
for the Gaussian curvature to be zero. This property can be used as a measure. For each 
resulting Bezier patch R(u,?;) = Il]=o 1(^)5 ^he Gaussian curvature, K, 

can be expressed in terms of ratios of bivariate Bernstein polynomials as [31] 



where S = R„ x Ry and S = |Ry x Ry|. If R(u,u) is a developable surface, then Kij must 
be zero for i = 0, . . . , 6n — 4 and j = 0, . . . , 2. However, for surfaces where F{x) 7 ^ 0, the 
Kij coefficients are not necessarily always precisely zero. 

From equation (5.37), we have 



K{u,v) 



(S • Ruu){S • Ryy) - (S . R„„)2 
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^denomiV-i 






(5.37) 




Knumj'^j l-^numC^i ^)| 



(5.38) 



^denom('^i'^) \7^d enomiV'1 '^)| 



\enom 



We can easily find the upper bound of |A"7ium('^? ^)| ^ 



6n-4 2 




2=0 j =0 
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6n— 4 2 
1=0 j =0 

6n-4 2 

< impn\Kij\ ^i,6n-4(w)^j,2('^) 

i=0 j=0 

= max|/fjj|. 
ij 

The denominator Kdenomi'^i'^) is always positive if the surface is regular, since x R^| 
is always positive for a regular surface. However, some of the coefficients Sij can be zero 
or negative. If there are non-positive Sij we can subdivide the surface into smaller patches 
until all Stj become positive. This is always possible by virtue of the convergence of the 
control polyhedron to the surface with repeated subdivision [22]. Then if all Sij > 0 we 
have the inequality 



\^denom{'^j'^)\ — 



> 



8n-4 4 

y^ 4(u)Rj^4(t^) 

t=0 j=0 

8n-4 4 

min|5,j| V V Bj,8n-4(«)%4(v) 
i=0 j=0 

min 5jj. 
u 



Therefore we can rewrite equation (5.38) as 



\K{u,v)\< 



maxjj \Kjj\ 

miuijSij 



(5.39) 



If minjj Sij is very small, the bound of the Gaussian curvature can be extremely large. 
In such case we further subdivide the patch and re-evaluate the bound until such bound 
converges to a stable value. 

Alternatively, the curvature may be analyzed at discrete locations on the surface, which 
provides a rough estimate but may not find the actual bound of the maximum curvature. 

We use the Gaussian curvature bound to give a physical measure of the developability of 
the surface. Throughout this thesis we say the surface is exactly developable if \K\ < 10“^^. 
For example, the curvature of the surface shown in Figure 5-4, as determined using ten sub- 
patches per Bezier patch in equation (5.39), is zero, so the surface is exactly developable. 



5.3.4 Developable Surfaces Constrained by Planes 

In many engineering applications, one or both of the directrices are restricted to planar 
curves which are not necessarily parallel. In this situation, we can force the control points 
hi and thus the directrix vb{u) to lie in a plane. This plane is determined by the two end 
points, bo and hmn for a degree (n, 1) B-spline patch with m Bezier segments, and a point 
g chosen by the user. The equation of the plane is given by 

ax + by-\-cz-\-d = 0 (5.40) 
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where, for a degree (n, 1) B~spline patch with m Bezier segments, 



a = {f>0 - g^Wmn - 9^) - iK - - 9^) 

b = {bl-9l{bl,n-9^)-ib^o-9^Wmn-9') 

C = ib^o-9^}{bln-9n-{bl-9^){bl,,,-9^) 

d = —ag^ — bg^ — cg^ . 



For each control point to lie on the plane, its rr, y and 2 components must satisfy 
equation (5.40). These lineax constraints can be used to reduce the number of unknowns in 
the equation using the generalized elimination method described in Section 5.3.1. 

However, if the (mn — 1) planar constraints, the (6m— 6) and continuity conditions 
and the two linear developability conditions £"0 and £m(3n-2) used to reduce the 

number of variables, the problem may become so overconstrained that is is difficult or 
impossible to find a feasible solution. For example, a cubic B-spline design curve with 
two Bezier segments has 15 unknowns (bi, b 2 , ba, b 4 , bs) and 13 linear equations, leaving 
just 2 variables for the entire problem. In this case, we use a subset of the total available 
linear equations as linear constraints and add the remaining equations to the objective 
function after squaring them. In selecting the linear equations to be used as constraints, 
it is important to remember that the linear constraints may be satisfied exactly while 
the objective function equations may not be. For example, if continuity is a must, 
equations (5.34) and (5.35) are used as linear equality constraints and the equations for 
planar constraints are included in the objective function. 





(a) Developable Surface with Constraining Plane (b) Final Degree (3,1) Developable Surface 



Figure 5-5: Planar Constraint on Developable B-Spline Surface Design 

Figure 5-5 shows a developable B-spline surface that was designed with the second 
directrix constrained to a plane. The same design curve and endpoints were used as those 
used for Figure 5-4, except the second directrix was constrained to lie in the z = 3 plane. 
Note that since the second directrix was constrained to a plane, the trivial solution was not 
an option, and simple bounds were not required. The continuity conditions were used as 
linear constraints. 

5.3.5 Degenerate Developable Surfaces 

Sometimes topologically triangular patches are required to produce a desired surface design. 
A three sided patch parameterized over a rectangular domain can be constructed by allowing 
one of the boundary generators to be zero length, say R(0, v). Faux and Pratt [11] showed 
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that the surface normal at this point may be evaluated using a limiting process based on a 
Taylor series expansion 



lim N(u, t>) 



Ku{0,v) X R^^(0,t;) 



For a degree (n-1) developable Bezier patch, equation (5.41) reduces to 



(5.41) 



N{0,v) 



(Pi,o ~ Pq,o) X (Pi,i ~ Pi,o) 
|(Pi,o-Po,o)x(Pi,i-Pi,o)|* 



Therefore, the unit normal vector at the degenerate point always exists as long as the three 
control points, Po,Oj Pi,o and Pi,i, do not coincide with each other, and it is orthogonal to 
the plane formed by the three control points as shown in Figure 5-6. 

The Gaussian curvature at the degenerate point cannot be calculated using the method 
described in Section 5.3.3 because the denominator, |Ru x vanishes. Wolter and 

Tuohy [43] devised a method for determining Gaussian curvature at a degenerate point using 
the second order partial derivatives of a local height function representation of the surface 
assuming the surface has a tangent plane and well defined curvatures at the degenerate 
point. For this description, assume that the degenerate point occurs at R(0,t>). The 
coordinate system is shifted such that the origin coincides with the degenerate point, the 
z axis is aligned with the normal to the surface and the x axis is aligned with the tangent 
to the surface in the non-degenerate parameter direction; in this case, R^(0, 0). The y 
axis is orthogonal to the others to form a right-handed system. Three linearly independent 
isoparameter curves in u are selected on the surface such that a„(u) = R(u, t;n )5 ^ = 
1, . . . , 3. The normal curvature in the direction 



^ aR(0,t;n) 
du 



^o;n(0) 

du 



{^n,yn) 



at the point R(0, 0) is given as 



n-k„ = no!"(0) 

= 9xx{0, 0)xl + 2gxy{0, 0)Xnyn + Oyyyl 

which can be represented in matrix form as 

[L]a = b 



(5.42) 



where 



^ “ (Qxxj 9xy^ 9yy) 

b = (n-a"(0),n-a2(0),n-a3(0))^ 



and 



[L] = 



xj 2xiyi y\ 

A 2x2 y2 yl 

A 2 x 3 y 3 y| 



If the fundamental forms are converted to the new coordinate system, the Gaussian curva- 
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ture can be represented as 



^ — 9xx9yy 9xy’ (5.43) 

The Gaussian curvature at a degenerate point can thus be determined by solving (5.42) for 
a and substituting the resulting values for Qxx^ 9xy and Qyy into (5.43). 

To estimate the maximum absolute value of the Gaussian curvature bound of a degen- 
erate surface with the degenerate point at R(0,0), the surface is first subdivided at u = e, 
where e is a small positive number so that the resulting patch with e < u < 1 is a nonde- 
generate patch. The Gaussian curvature is calculated at w = 0 using equation (5.43), and 
the bpund for the Gaussian curvature for the nondegenerate patch is evaluated using the 
technique described in Section 5.3.3. 



foi 




Figure 5-6: Triangular Degenerate Cubic Bezier Patch 

Figure 5-6 shows a degenerate cubic-linear Bezier patch with non-planar directrices. The 
design curve is made up of the control points (0, 0, 0), (1.8, 3, -0.7), (3.3, -2, 1.5) and (4, 
0.5, 0). The selected end points are (0, 0, 0) and (4, 0.5, 2). The Gaussian curvature at the 
degenerate point was found to be 8.91 x 10”^"^. The surface was split at the parameter value 
u — 0.0001 resulting in a maximum bound for the Gaussian curvature of zero. Therefore, 
the surface is exactly developable. Note that there are only three points which define the 
end rulings and they always lie on the same plane. Therefore, there is no skew between the 
two end rulings and the surface often results in an exact developable surface. 



5.4 Development of a Developable Surface onto a Plane 

In the manufacture of developable surfaces, it is necessary to calculate the plane development 
of these surfaces. This determines cut information for the flat piece of metal so it can be 
cut to the proper dimensions before it is rolled. 

This formulation follows the method of Faux and Pratt [11, p.274], using the fact that a 
developable surface may be mapped isometrically onto a plane. The Frenet-Serret formulae 
state that 



t' = «n 



(5.44) 
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where t is the unit tangent vector to a curve, n is the unit normal vector to a curve and 
K is the curvature. For a planar curve in the (X, Y) plane, we can define the unit normal 
vector as 



n = Cz X t 



(5.45) 



where = (0,0, 1). Substituting (5.45) into (5.44) yields 

(fX dY ^ 

<PY dX „ 



(5.46) 

(5.47) 



where (X, Y) denote the 2D coordinates on the developed plane (X,Y). If we rewrite equa- 
tion (5.46) in terms of the parameter u, we obtain 



(fX ^ 

ds^ ds 



ds \ ds J 



' ds 



__ d / dX du \ dY du 

ds \ du ds J du ds 



dX f(fu 

du"^ \ds ) du i 



^ du ds 



If {du/ds)*^ ^ 0, this reduces to 

d^X 

du^ 



d^Y 

du^ 



+ 



iv) 


dX 


K 


dY 




du 


(S) 


du 


represented 


as 




(&) 


dY 


n 


dX 


ay 


du 


(s) 


du 



= 0 . 



(5.48) 



= 0 . 



(5.49) 



The development is based on the fact that curves on isometric surfaces have the same 
geodesic curvatures. It can be seen from Figure 2-4 that for a planar curve, since the angle 
between the surface normal and the normal to a curve is 90 degrees, the curvature vector 
equals the geodesic curvature vector. Therefore, k in equations (5.48) and (5.49) can be 
replaced by the geodesic curvature of the curve on the developable surface. If we choose 
the curve on the developable surface to be an isoparametric curve in terms of u, R(u,Un), 
we can replace du/ds and d^u/ds^ in equations (5.48) and (5.49) by 



du 

ds 

d^u 

ds"^ 



1 



and 



|R^(u, Vfi)\ 

Rtz(u,Un) • Riiix(u,Vn) 
Ru(u, Vn) * Rit(u, Vn) 
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Thus we have 



du 

du 

dp 

du 

dq 

du 



P 



Q 






(5.50) 



The surface will be developed onto the plane oriented such that the point R(0,0) is 
located at the origin of the {X^Y) coordinate system in the plane and the vector Ru(0,0) 
coincides with the X axis, as shown in Figure 5-7. We integrate the system (5.50) along 
the directrix that corresponds to i; = 0 from u = 0 to w = 1, using, for example, a Runge 
Kutta scheme. Since isometric maps are conformal, the angle between the directrix and the 
generator at (0, 0) is the same in both representations and can be found by 

. Rt.(0,0) R.(o,o) 

|Ru (0,0)1 ■ |Rv(0,0)|- 

Straight lines on the surface map to straight lines on the plane, so the generator length 
on the surface is measured and transferred to the plane, giving the initial point for the de- 
velopment of the second directrix. Then the system of ordinary differential equations (5.50) 
is integrated using a fourth order Runge-Kutta method [38] along the directrix that cor- 
responds to V = 1 from u = 0 to u = 1. Finally, the endpoints of the two directrices are 
connected. 

To verify that the development of the surface is accurate, one can check the results 
by computing the angle between the directrix corresponding to v = 0 and the generator 
corresponding to u = 1, computing the length of that generator, and determining the 
location of the point R(l,l) in the plane. Comparing this point to the point calculated 
above determines the accuracy of the development. The development will be more accurate 
with increased steps in the Runge-Kutta scheme. 

The surface that was designed in Section 5.3.2 is developed onto a plane here. This 
surface is shown in Figure 5-7(a). It was developed onto a plane using 31 steps in a fourth 
order Runge-Kutta routine as shown in Figure 5-7(b). The accuracy for the check at R(l, 1) 
was 10“^. 



5*5 Examples 

5.5.1 Example of a Developable B-Spline Strip Surface 

This example is a windshield of an automobile as depicted in Figure 5-8. Since automobile 
windshields are generally sandwiched between two parallel planes, we can use the technique 
introduced in Section 5.2. The design curve, r^(u), is given in the xy-plane by a cubic 
B-spline curve with control points (0, 2.1, 0), (0.5, 0.9, 0), (1.5, 0.3, 0), (4.5, 0, 0), (7.5, 0.3, 
0), (8.5, 0.9, 0), (9, 2.1, 0) and a knot vector (0, 0, 0, 0, 0.3, 0.5, 0.7, 1, 1, 1, 1). The two 
end boundary conditions are chosen to be bo = (0.2, 3.5, 1.5) and be = (8.8, 3.5, 1.5). The 
design curve is first split into four cubic Bezier curves and the control points of r^(u) are 
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R(0,0). 




(a) Developable B-Spline Surface 



(b) Plane Development of Surface 



Figure 5-7: Development of B-Spline Developable Surface 



evaluated. The resulting B-spline curve, ra(w), has internal control points (0.56, 2.64, 1.5), 
(0.96, 2.18, 1.5), (1.84, 1.81, 1.5), (3.27, 1.62, 1.5), (4.5, 1.53, 1.5), (5.73, 1.62, 1.5), (7.16, 
1.81, 1.5), (8.04, 2.18. 1.5) and (8.44, 2.64, 1.5), with knot vector (0 0 0 0 0 .3 .3 .5 .5 .7 
.711111). Note that the resulting surface is not cylindrical. After degree elevating the 
design curve by one and removing the knots from T 5 (u), we obtain the degree (4-1) B-spline 
developable surface shown in Figure 5-8. The windshield surface is then developed onto the 
plane as shown in Figure 5-8 (c). 

5.5.2 Example of a Developable B-Spline Surface with 3D Directrices 

The stack of a ship, depicted in Figure 5-9, was constructed using a B-spline design curve 
with the resulting curve constrained to a plane not parallel to the plane of the design curve. 
Therefore, the method described in Section 5.3 was required. The stack was split along the 
plane of symmetry (the xz plane) from forward to aft. The simple shape allowed definition 
of the design curve with a single B-spline curve of two Bezier segments for each half stack. 
A cubic design curve was selected that defined the base of the stack, with control points (0, 
0, 0), (0, 10, 0), (15, 10, 0), (35, 3, 0) and (35, 0, 0) and knot vector (0, 0, 0, 0, .5, 1, 1, 1, 
1). The two endpoints were selected to be bo = (6, 0, 31) and be = (24, 0, 29). The second 
directrix was constrained to lie on a plane, with the third point to define the plane being 
(24, 10, 29). The optimization process was conducted using equations (5.34) and (5.35) to 
reduce the number of variables, and no simple bounds were employed. The resulting surface 
was degree (3,1). The maximum Gaussian curvature determined using equation (5.39) with 
10 Bezier subpatches was 2.19 x 10“^. The surface was then reflected with respect to the 
xz plane, resulting in a symmetric body. One of the surfaces surfaces was developed onto 
a plane, as shown in Figure 5-9 (c). 
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(a) Developable Surface with Control Points 





Figure 5-8: Automobile Windshield, a Developable Degree (4-1) B-spline Strip Surface 
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Figure 5-9: Ship Stack, A Developable Degree (3,1) B-Spline Surface with \K\ < 2.19 x 10 ^ 
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Chapter 6 

Geodesics on Developable Surfaces 



6.1 Introduction 

Geodesics axe useful in ship design and shipbuilding for determining the layout of seams and 
butts in a ship hull, as described by Munchmeyer and Haw [33]. This chapter shows that a 
two point boundary value problem (BVP) for computing geodesics on a developable surface 
can always be reduced to an initial value problem (IVP), with a corresponding savings in 
time for the solution process. Most problems that arise in applications of geodesics are 
BVPs, which are much more difficult to solve than IVPs. 



6.2 Formulation 



Recall from Section 2.3 that the governing equations of a geodesic are given by a set of 
coupled second order ordinary differential equations (2.29) and (2.30) 



— + ri 



(fiv 

ds^ 



73 + ^^11 




+ 2T 



1 dudv 1 




9 du dv 9 



22 




= 0 
= 0 



where u and v are related by the first fundamental form ds^ = Edv? + 2Fdudv + Gdv^. 
These two second order differential equations can be rewritten as a system of four first order 
differential equations 



du 

ds 

dv 

ds 

dp 

ds 

dq 

ds 



P 

Q 

-T\,p'^-2T\2Pq-VW 

-ri,p^-2r\^pq-rl^q^. 



( 6 . 1 ) 

( 6 . 2 ) 

(6.3) 

(6.4) 



We can solve the system as an initial- value problem (IVP), where all four boundary 
conditions are given at one point, or ais a boundary- value problem (BVP), where the four 
boundary conditions are specified at two distinct points. Most of the problems that arise 
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in applications of geodesics are BVPs, which are much more difficult to solve than IVPs. 

The solution of an IVP is unique; however for a BVP it is possible that the differential 
equations will have many solutions or even no solution [25], General methods for the 
solutions of two-point BVPs can be found in [12, 25]. There are two commonly used 
approaches to the numerical solution of a BVP, namely the shooting method and the finite 
difference method (relaxation method). The shooting method was originally used to adjust 
the settings of artillery equipment to hit a target. To implement the method, the unknown 
boundary conditions are guessed at the initial point and an IVP method is implemented. 
Based on the resulting values at the end point, the initial guess is adjusted until the “target” 
is “hit”. The second approach is based on a finite diflFerence approximation to ^ where 
y = {u,v,p,q) on a mesh of points in the interval [A,B]. This method starts with an 
initial guess and improves the solution iteratively. Maekawa [30] solves the geodesic BVP 
on free-form parametric surfaces by the finite diflFerence method. In general, the shooting 
method is very sensitive to the initial guess at point A, and is much less stable than the 
finite difference method [30]. 

For a B-spline developable surface patch there is a single unique solution to the system. 
Since a geodesic on a developable maps to a straight line on a developed plane, there is 
only one solution to the system on a developable surface. Here we exclude periodic surfaces 
such as cylinders where there can be more than one solution, as shown in Figure 6-1. 




Figure 6-1: Geodesics on a Closed Surface 

In this section it is shown that all two point BVPs for solving geodesics on developables 
can be reduced to IVPs using the following properties described in Section 4.2: 

• A developable surface can be mapped isometrically onto a plane. 

• Isometric surfaces have the same Gaussian curvature at corresponding points. 

• Corresponding curves on these surfaces have the same geodesic curvature at corre- 
sponding points. 

• Every isometric mapping is conformal; i.e. the angle of intersection of every arbi- 
trary pair of intersecting arcs on a developable surface is the same as that of the 
corresponding inverse image in the plane at the corresponding points. 
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• A geodesic on a developable surface maps to a straight line in the plane. 

The basic procedure is to map the two desired points on the developable surface to a plane, 
draw the straight line between them and determine the angle between the generator and 
the geodesic line at one of the end points. The angle can be used to determine du/ds and 
dv/ds. Thus, all the information required for an IVP is available. 





(a) Surface Developed onto Plane 



(b) uv Parametric Space 




Figure 6-2: Geodesic on a Degree (3,1) Developable Surface 

Given two points A and B on the developable surface R(u, u) as shown in Figure 6-2(c), 
the corresponding points {Xa^ Ya) and {Xb^ Yb) in the developed plane are required. To find 
the points in the plane, (Xq, To) is set as the (0, 0) point in the plane corresponding to R(0, 0) 
on the surface. Using the method presented in Section 5.4, the directrix corresponding to 
u = 0 is developed into the plane to determine the point C = (u>v,0), shown in Figure 6- 
2(a). The angle between the directrix and the generator at (u^,0) is the same in both 
representations and can be found by 



cos d = 7—7 -77 r- 

\Ru{uaM i^i 



where is an isoparametric line R(u>i, u) which is a straight line on the surface. Therefore 
it is a geodesic and will be developed into the plane as a straight line. The distance \C^\ 
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is given by 



1^1 = \J {XA - + (j/A - 2/c)^ + {ZA - Zc)'^. 

The point A on the plane is found using C, \C^\ and 6. The point B on the plane is found 
by following the same procedure, and the points are connected as shown in Figure 6-2(a). 

The angle u between cX and is given by 



COSCJ = 



\Ai\\cX\ 



This angle u is preserved between the isoparametric line t{ua,v) and the geodesic curve 
g(s) on the developable surface at point A. Thus we have 



cos a; = 



|Rn,||g'(^)| 



where the tangent vector to the geodesic is given by 

//X ^ du ^ dv _ 

g (s) — -h 

The angles cj and 9 are shown in Figures 6-2(a) and (c). Multiplying (6.5) by Ry yields 



• g'(5) = R^ • R^— -f Ky • Ry— = cos(u;)|Ry||g'(s)| 
as as 

which (since |g'(5)| = 1) can be reduced to 

^du ^dv , , / 7 - 
F— + G— = cos(u;)vG 
ds ds 

where F = Ry • Ry and G = Ry • Ry (coefficients of the first fundamental form). Thus, 



dv _ cos{u) F / du\ 
ds y/G G \ds) ‘ 



From the first fundamental form. 



g'(5).g'(s) = S 




^^dudv ^ 
+ 2F— — + G 
ds ds 




= 1 . 



Plugging (6.6) into (6.7) and solving for du/ds yields 



and thus (6.6) reduces to 



du _ ^ lsin^{u)G 

EG- F^ 



dv _ cos(ct;) ^1 F 
^ ~ y/G 7^ y/E G - ^ 



(6.6) 



(6.7) 
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Hence, we have all the initial conditions required to solve the IVP (equations (6.1) to (6.4)) 
for a geodesic. The solution to the IVP yields the uv parametric values for the geodesic that 
are graphed in Figure 6-2 (b). The corresponding three-dimensional coordinate values are 
shown in Figures 6-2(c) and (d). The developable degree (3,1) surface shown in Figure 6-2 
has the control points and knot vectors shown in Table 6.1. The geodesic runs from 

{ua.va) = (0.1, 0.3) to [ub^vb) = (0.9, 0.8). 



Po,0 

Pl.O 

P2,0 

Ps.o 

P4,0 


(-2, 1, 0) 
(0, -1, -1) 
(2, 5, -1.4) 
(4, 1, -1) 
(6, 5, 0) 


Po,i 

Pi,o 

P2,0 

Pa.o 

P4,0 


( 0, 2, -3) 
(1, 1, -3.5) 
(2, 4, -3.7) 
(3, 2, -3.5) 
(4, 4, -3) 


Tu 


0, 0, 0, 0, 
.5, 1,1, 1,1 




0, 0, 1, 1 



Table 6.1: Control Points for a Degree (3,1) Developable Surface 



6.3 Examples 

In this section we conduct numerical experiments to compare the performance of the IVP to 
the finite difierence based BVP on developable surfaces. For a fair comparison, we determine 
the minimum number of steps in the BVP to reliably meet the given tolerance e = 10“^ 
for Newton’s method, see [31]. We then find the minimum number of steps for the IVP 
to achieve the same order of accuracy at the endpoint B. The examples were run on a 180 
MHz Silicon Graphics workstation. 

The example shown in Figure 6-2 runs using an IVP solution with 66 steps in 0.25 
seconds, as compared to 0.60 seconds for a BVP solution with 135 steps and 7 iterations. 
This is over twice as fast, with a very simple surface and a geodesic that is very close to a 
straight line in uv space. 

A more complex degree (4,1) developable surface with five Bezier patches also runs over 
twice as fast: 0.88 seconds for the IVP solution with 194 steps instead of 2.28 seconds for 
the BVP solution with 294 steps and 11 iterations. The control points and knot vectors for 
this example are displayed in Table 6.2. The geodesic runs from {ua^va) = (0.1, 0.9) to 
{ub^vb) = (0.9, 0.1). The surface and its geodesic are displayed in Figure 6-3. 

The final example is a handkerchief-like surface. This degree (4,1) developable surface 
with four Bezier patches runs in 0.50 seconds for the IVP and 2.66 seconds for the BVP, 
which is over 5 times faster. The IVP required 176 steps, while the BVP required 320 
steps and 12 iterations. The control points and knot vectors for this example are shown in 
Table 6.3. The geodesic runs from (ua^va) = (0,0) to {ub,vb) = (ljl)» The surface and 
its geodesic are displayed in Figure 6-4. 
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(a) Surface Developed onto Plane 



(b) uv Parametric Space 




(c) Geodesic on 3D Surface 

Figure 6-3: Geodesic on a Degree (4,1) Developable Surface 



Po,o 


(0.00, 3.00, 0.00) 


Po,i 


(0.00, 3.00, -5.00) 


Pl,0 


(0.75, 1.50, 0.00) 


Pi,i 


(2.42,-1.84, -5.00) 


P2,0 


(1.25, 1.50, 0.00) 


P2,l 


(3.88,-1.53, -5.00) 


Pa.o 


(1.96, 2.58, 0.00) 


P3,l 


(5.85, 1.34, -5.00) 


P4,0 


(2.50, 2.00, 0.00) 


P4,l 


(7.12,-2.28, -5.00) 


Ps.O 


(3.00, 1.17, 0.00) 


Ps.i 


(8.21,-1.97, -5.00) 


Pe.o 


(3.50, 1.00, 0.00) 


Pe.i 


(9.09,-2.20, -5.00) 


P7,0 


(4.00, 1.08, 0.00) 


P7,l 


(9.86,-2.03, -5.00) 


Pg.o 


(4.50, 1.50, 0.00) 


Ps,i 


(1.04,-1.53, -5.00) 


Pg.o 


(5.04, 2.29, 0.00) 


P9,l 


(1.09,-6.81, -5.00) 


Pio,o 


(5.75, 4.25, 0.00) 


PlO,l 


(1.13, 4.04, -5.00) 


Pll,0 


(6.75, 4.50, 0.00) 


Pll.l 


(1.18, 1.05, -5.00) 


Pl2,0 


(9.00, 3.00, 0.00) 


Pl2,l 


(1.20,-1.11, -5.00) 




0, 0, 0, 0, 0, .2, .2, 
.4, .4, .6, .6, .8, .8, 


Ty 


0, 0, 1, 1 




1, 1, 1, 1, 1 







Table 6.2: Control Points for a Degree (4,1) Developable Surface 
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(a) Surface Developed onto Plane 




(b) uv Parametric Space 




Figure 6-4: Geodesic on a Degree (4,1) Handkerchief-like Developable Surface 



Po,o 


(0.50, 3.00, 0.00) 


Po,i 


(0.50, 3.00, -5.00) 


Pl.o 


(0.88, 1.50, 0.00) 


Pl.o 


(0.88, 1.50, -5.00) 


P2,0 


(1.25, 2.00, 0.00) 


P2,0 


(1.25, 2.00, -5.00) 


P3.O 


(1.96, 4.17, 0.00) 


Pa.o 


(1.96, 4.17, -5.00) 


P4,0 


(2.50, 3.00, 0.00) 


P4,0 


(2.50, 3.00, -5.00) 


Ps.o 


(3.00, 1.67, 0.00) 


Ps.O 


(3.00, 1.67, -5.00) 


Pe,o 


(3.50, 3.00, 0.00) 


Pe.o 


(3.50, 3.00, -5.00) 


P7,0 


(4.04, 4.17, 0.00) 


Pr.o 


(4.04, 4.17, -5.00) 


Ps,0 


(4.75, 2.00, 0.00) 


Ps,0 


(4.75, 2.00, -5.00) 


Pg,o 


(5.13, 1.50, 0.00) 


Pg.o 


(5.13, 1.50, -5.00) 


Pio,o 


(5.50, 3.00, 0.00) 


Pio,o 


(5.50, 3.00, -5.00) 


T„ 


0, 0, 0, 0, 0, 


Tv 


0, 0, 1, 1 




.25, .25, .5, .5, .75, .75, 








1, 1, 1, 1, 1 







Table 6.3: Control Points for a Degree (4,1) Handkerchief-like Developable Surface 
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Chapter 7 

Lines of Curvature on Developable 
Surfaces 



7.1 Introduction 

This chapter discusses the behavior of the lines of curvature on a developable surface near 
the inflection line. A developable surface does not possess isolated generic flat points but 
rather may contain a line of non-generic flat points along the generator. This property can 
be used to detect an inflection line where the (only) non-zero principal curvature whose 
principal direction is perpendicular to the generator changes its sign. Munchmeyer [33] 
indicates that a developable surface can be shaped purely by rolling and should be fed into 
the roller so that the direction of the generator is parallel to the rolls. However, when the 
non-zero principal curvature changes sign, the plate can no longer be fed into the roller 
in the same direction, since the bending direction changes at the inflection line. Locating 
the inflection line prior to fabrication allows the manufacturer to either cut the plate along 
the inflection line or reverse the direction the plate is fed through the rollers at that line. 
Surface inflection was studied by Hoitsma [21], who showed that a surface has an inflection 
at a point P if and only if its mean curvature changes sign in the neighborhood of P. This 
chapter further extends this result. 



7.2 Inflection Line 



This section establishes the definition of an inflection line on a developable surface. First, 
some additional differential geometry properties are required. 

Recall that throughout this thesis, it is assumed that the u = const isoparametric line 
corresponds to the generator of the developable surface or, in other words, the straight line 
ruling is in the v direction. With this assumption, Hyy = 0; consequently, from (2.13), 
N = Rvt, • N = 0. From (2.23), since Gaussian curvature of a developable surface is zero. 



LN-M^ _ -M^ 

EG -F^ ~ EG- F2 



(7.1) 



and we have M = 0. Therefore, mean curvature as expressed by equation (2.24) reduces to 



H = 



2FM -EN-GL 
2{EG - F2) 



-GL 

2{EG-F^y 



(7.2) 
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Recall from (4.22) that the nonzero principal curvature is given by = 2H, When L = 0, 
H and k* become zero; otherwise, k* = 2H 7^ 0. 

Next, we will show that the u = const parametric straight lines become the lines of zero 
curvature. This can be seen from the fact that 

• the u = const paxametric straight lines have zero normal curvature, and 

• no other direction has zero normal curvature. 

The second fact comes from Euler’s theorem (2.25), 

Kn = Kicos^(a) + «2sin^(a). 

When K 2 = 0, cos^(a), which becomes zero only when a = | or corresponding 

to the direction of Similarly, when = 0, = 0 only when a = 0 or tt, corresponding 

to the direction of k\. 

Theorem 7*2. 1 A developable surface does not possess generic isolated flat points ^ but 
rather may contain a line of non~generic flat points along a generator. 

Proof: From equation (7.2), L vanishes at a flat point {u/^Vf) where K = H = 0. There- 
fore, from equation (2.14) we have 

L{uf,vf) = -Ru{uf,Vf) * Nu{uf,Vf) = 0. (7.3) 

From equation (7.1), M = 0 on a developable surface. Therefore, from (2.15) we have 

M{uf,vf) = ~R^{uf,Vf) • Nu{uf,vf) = 0. (7.4) 

Since N is a unit vector, we also have 

N(uf,Vf) • Nu{uf,Vf) = 0. (7.5) 

If Nti is not zero, then from (7.3), (7.4) and (7.5) must be perpendicular to Ru, 
and N. This is impossible because N is perpendicular to both R^ and R^;, and R^ is not 
parallel to R^. Thus, Nu(u/,u/) must equal zero. 

For developable surfaces, the unit normal vector N is constant along a generator. There- 
fore, the rate of change of the unit normal vector in the u direction must also be constant 
along a generator. This leads us to the fact that Nu is not only zero at {uf,Vf) but also 
zero along the u = Uf isoparametric line. Therefore, for a given u = Uf^ equation (7.3) 
becomes 



L{uf,v) = -R^^u(^i/,?;) • N{uf,v) = 0 (7.6) 

for 0 < 1; < 1. Consequently, the entire generator consists of a line of flat points. I 

For a developable surface, the inflection line is a generator which consists of a line of 
flat points and coincides with a change in sign of the nonzero principal curvature. In other 
words, the nonzero principal curvature has a different sign on either side of an inflection 
line. 



* A developable surface can not possess sphericed umbilics since one of the principzd curvatures is always 
zero. 
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The inflection line can be detected by finding u = uj such that L{uj^Vn) = 0 where 
is an arbitrary constant between 0 and 1. L{u^Vn) = 0 can be written as 



Ruit (u, Vn) 



Vn) X q 

|R^(U,U„) X R^(u,Un)| 



(7.7) 



Since we are atssuming a regular surface such that |Rtx x R-vl 7^ 0, we need only set the 
numerator of equation (7.7) to zero. Thus, 



|Ruu('^»'^n) R^('^j'^n) R^('^7'^n)l — 0- 



For a polynomial surface with degree n in the u direction, this results in a univariate 
polynomial equation of degree (3n — 4) in u 



{Vu^v ^uyv)^uu (^v^u ^u^vJPuu i^uUv ^vyu)^uu — 0 . 



(7,8) 



If the surface is expressed in a piecewise polynomial form such as a B-spline representation, 
equation (7.8) must be applied to each polynomial segment separately. The univariate poly- 
nomial equation can be robustly and efficiently solved by the Interval Projected Polyhedron 
algorithm [24, 39]. Once we find the u value, we check if the nonzero principal curvature 
changes sign. If the sign does not change, the generator is not an inflection line. A simple 
example is given by R(u,u) = (u, u, u"^). Since 

G=l, Ru„ = (0,0,12 u 2) and N = 



we have 



, — 12u^ 

^ 

\/l6u^ + 1 

It is apparent that k* becomes zero when u = 0; however k* does not change sign when u 
moves from u < 0 to u > 0. Therefore, u = 0 is not an inflection line. 

At flat points the principal directions are indeterminate and the orthogonal net of lines 
of curvature may have singular properties. In the following we investigate the pattern of 
the lines of curvature near the line of non-generic flat points. Locally any surface is the 
graph of a differentiable function [8]. Given a point p of a surface 5 we can choose the 
coordinate ajcis so that the origin O of the coordinate is p and the z axis is directed along 
the positive normal of S at p; thus, the xy plane coincides with the tangent plane of 5 at 
p. It follows that a neighborhood of p in 5 can be represented in the form z = h{x,y) with 

/i(0,0) =/ix(0,0) =/iy(0,0) =0. 

By taking into consideration that /i(0, 0) = /ix(0, 0) = /iy(0, 0) = 0, the Taylor expansion 
of the z component of the surface becomes 

^(^1 y) = 0 ) + 2xyhxy{0, 0 ) + y'^hyyiO, 0 )] 

+ ^[3:^/ixxx(0,0) + 3a:^j//ixxy(0,0) + 3xy‘^hxyy{0,0) + y^hyyy{0,0)] 

+i?(x,y)(l2:,y|^) (7.9) 

where R{x, y) is a remainder term with limx_+o,y->o y) = 0 and |rr, y| = \/x'^ + y^. 
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This local approximation will now be applied to developable surfaces. 

Lemma 7.2.1 A developable surface is, in general, locally a parabolic cylinder and becomes 
a cubic cylinder at the inflection line. 

Proof: Let us consider an orthogonal Cartesian reference frame attached to the surface at 
an arbitrary point (woj^o)* We choose unit vectors x N, and N as the directions 
of X, y and z axes such that the y axis coincides with the generator u = uq, the z axis 
coincides with the normal vector and the x axis is orthogonal to both axes. Equation (7.9) 
reduces to [32] 



y) = 0) + 0) + 3x^y/iiiy(0, 0)] + R{x, y){\x, (7.10) 

Using the inverse function theorem, /ixx(0,0), /ixxx(0?0) and /ixxy(0, 0) can be obtained as 
a function of u and v. At a the line of inflection, /ixx(0,0) and /ixxi/(0, 0) become zero [32], 
and equation (7.10) reduces to 

h{x, y) = i(N • x^ + R{x,y){\x,y\^) (7.11) 

where the coefficient of x^ is evaluated at (uo,t>o)- From equation (7.10) it is apparent that 
for small x, the quadratic term dominates. From equation (7.11) it is apparent that at an 
inflection line, the surface is locally a cubic cylinder. 1 

Theorem 7.2.2 There is only one line of curvature that passes through each fiat point on 
an inflection line, and that line of curvature is orthogonal to the direction of the generator. 



Proof: By Lemma 7.2.1 the developable surface is expressed locally as a cubic cylinder. 
If we rewrite the cubic term of equation (7.11) in terms of polar coordinates by substituting 
X = r cos 0 we obtain 



h{e) = j{N.r,,u) cos^e. ( 7 . 12 ) 

Since h{9 + tt) = —h{9), h{9) is an antisymmetric function of 9, The roots of ^ = 0 will 
give the angles where local maxima and minima of h{9) may occur around the flat point. 
The equation can be restricted to the range 0 < 9 < 27t without loss of generality. The 
roots are easily computed as 0 = 0, |, tt and Only 9 = 0 and 9 = ir (which coincide 
with the local x axis) give extrema, since 

d9^ d9^ 

and thus 0 = f and 9 = ^ (which coincide with the local y axis) provide neither a maximum 
nor a minimum. In other words, the inflection line is not a line of curvature. Consequently, 
there is only one line of curvature (minimum curvature where 9 = 0 and maximum curvature 
where 9 = tt) that passes through the flat point. I 
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Recall from Section 2.3 that the principal curvatures can be found by solving the two 
simultaneous equations (2.19) and (2.20), 

(L + KE)du + (M + kF)(1v = 0 
(A/ + KF^du -h (iV “H KG^dv = 0, 

for K. Since M = A/' = 0, equations (2.19) and (2.20) reduce to 

{L KE)du + KFdv = 0 and (713) 

KFdu-\~KGdv = 0. (7.14) 

We can trace the lines of curvature which pass through the flat points of an inflection line 
by providing the starting points and shifting outwards in the directions 0 and tt from the 
flat points or, equivalently, along the positive and negative local x atxis. The tracing can 
be carried out by integrating the initial value problem for the coupled ordinary differential 
equations (7.13) and (7.14). 

Generic umbilics (including flat points) are stable with respect to small perturbations of 
the function representing the surface, while non-generic ones are unstable [32]. In generic 
cases umbilics are isolated [20]; thus the inflection line, which consists of line of flat points, is 
non-generic and therefore unstable. In the following we give a couple of numerical examples 
that demonstrate the instability of the line of flat points along the inflection line with respect 
to perturbations. 

The example surface is a degree (3-1) integral Bezier patch which is constructed by the 
method described in Section 5.3, see Figure 7-1. The control points are shown in Table 7.1. 
The surface has an inflection line at = 0.5754, which has been computed by solving the 



Po ,0 


(0.0000, 0.0000, 0.0000) 


Po.i 


(0.5000, 0.0000, 2.0000) 


Pl,0 


(1.8000, 3.0000, 0.0000) 


Pi,o 


(1.8950, 2.3250, 2.0000) 


P2,0 


(3.3000,-2.0000, 1.5000) 


P2,0 


(3.0575,-1.5500, 3.1625) 


Pa.o 


(4.0000, 0.0000, 0.0000) 


Ps.o 


(3.6000, 0.0000, 2.0000) 



Table 7.1: Control Points for a Degree (3,1) Developable Surface 



degree 5 univariate polynomial equation (7.8). The inflection line is shown in Figure 7-1 
with a dash dotted line. This surface has a net of lines of curvature which is shown in 
Figure 7-2(a). Solid lines represent the lines of maximum principal curvature, while dotted 
lines represent the lines of minimum principal curvature. The inflection line is depicted with 
a dash dotted line. Figure 7-2(b) shows a magnification near the inflection line. We can 
observe that there is only one line of curvature that passes through a flat point orthogonal 
to the inflection line. 

We gradually perturb the control points of the surface and observe the behavior of the 
lines of curvature which pass through the inflection line. The control points are perturbed in 
the following manner. Since the example is a degree (3-1) patch, it has 8 control points. Each 
control point consists of three Cartesian coordinates rr, y, Z] hence there are 24 components 
to be perturbed. A random number which varies from —1 to 1 is used to determine the 
perturbation of each component. Let us denote the randomly chosen numbers for each 
control point as (e^-, ef^), 0 < i,j <3. We normalize the vector and add it to each 
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Figure 7-1: Developable surface and its control polyhedron with inflection line. 




Figure 7-2: (a) Lines of curvature of developable surface with inflection, (b) Magnification 
near inflection line. 
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Figure 7-3: (a) Lines of curvature on perturbed surface ^ = 0.02. (b) Magnification near 
u=0.57. 





Figure 7-4: (a) Lines of curvature on perturbed surface C = 0.08 (b) Magnification near 
u=0.57. 



87 



control point as follows: 



Pij — Pij + C 



[e 



u-i 






2 . eV.2 . 2 



(7.15) 



where C = 0.02. We gradually increase the perturbation by increasing ^ from 0.02 to 0.08 
in steps of 0.02. 

Figures 7-3 and 7-4 illustrate the behavior of the lines of curvature when the control 
points are perturbed. We can see from the figures that the entire inflection line, which 
consists of a line of flat points, disappears. Hence there is no singularity in the net of lines 
of curvature when a perturbation is induced. The nonzero principal curvatures on both 
sides of the former inflection line ^ meet at right angles near the former inflection line and 
make a very sharp change in direction (almost a right angle). 



^Once the control points are perturbed both principal curvatures may not be nonzero, but here we zu*e 
referring to the nonzero principal curvature before perturbation. 
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Chapter 8 

Engineering Example 



The final example is the design of a small boat using B-spline developable surfaces. The 
boat wats designed with a hard chine which separates the bottom and side of the boat and 
demarks a tangent discontinuity in the hull of the boat. As an initial stage of the design, the 
following three curves, which define the desired characteristics of the hull, were determined 
using a model of a planing boat: 

• sheer curve - the line that runs along the deck edge from bow to stern 

• baise curve - the line that runs along the keel of the boat (centerline from bow to 
stern) 

• chine curve - the line that runs from the bow to the lower corner of the stern 

Figure 8-1 shows the terminology used in this section. The coordinate system is located 
so that the origin is located on the baseline at the bow of the boat with the x axis pointing 
aft, the y axis pointing to starboard and the z axis pointing up. (Note that the baseline, 
a zero reference line, is different from the base curve defined above.) Due to the symmetry 
of the boat, we design only the starboard side and reflect those values across the plane of 
symmetry to obtain the port side. The cubic B-spline representation of the three curves is 
shown in Table 8.1. 




Figure 8-1: Planing Boat Terminology 



Design Using 3D Directrices 

Employing the design philosophy described in Chapter 5, the chine curve was designated 
cLS the design curve and the sheer and base curves as target curves, i.e. ra(u), to construct 
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Sheer Curve 


Chine Curve 


Base Curve 


Control 

Points 


(0.00, 0.00, 9.00) 
(6.86, 7.10, 8.22) 
(21.6, 8.93, 6.25) 
(36.9, 8.73, 5.86) 
(45.0, 7.65, 6.10) 


(1.40, 0.00, 5.30) 
(10.5, 7.53, 1.93) 
(25.7, 7.85, 1.28) 
(40.4, 7.46, 1.27) 
(44.1, 7.20, 1.70) 


[l.40, 0.00, 5.30) 
(2.26, 0.00, -0.21) 
(22.6, 0.00, -0.10) 
(36.3, 0.00, -0.10) 
(44.1, 0.00, 0.50) 


Knot 

Vectors 


0, 0, 0, 0, .5, 
1, 1, 1, 1 


0, 0, 0, 0, .75, 
1, 1, 1, 1 


0, 0, 0, 0, .5, 
1, 1, 1, 1 



Table 8.1: Target Curve Data 



the side and bottom respectively. Initially, the endpoints of the sheer and baise curves were 
used to construct the side and bottom of a single surface each; however, this resulted in a 
flat side and bottom. To remedy this situation, the chine was split into two B-spline curves 
at a parameter value of u = 0.5, and one additional knot wcls inserted at the midpoint of 
the forward B-spline curve. The resulting curves are shown in Table 8.2. Thus, the side 
and bottom are constructed of two surfaces each, allowing more flexibility in the design, 
but losing the continuity along the surface. The surfaces are, however, at least 
continuous between the forward and aft portions. Since the design curve wcLS split into two 
curves that maintain continuity between them, the tangent at the connecting point is 
continuous. Since both patches are developable, both surfaces share the same tangent plane 
at the connecting generator. Hence, the patches are continuous. 

Note that the forward and aft surfaces join along a generator which is a straight line. 
The planar transom (stern) joins the aft surfaces along generators as well. 

The endpoints for the side and bottom surfaces were selected so as to force the resulting 
directrix as close to the target curves as possible while maintaining a majcimum Gaussian 
curvature on the order of 10“^. The data for each developable surface on the starboard side 
of the boat is listed in Table 8.2. Measurements are in feet. The port side of the boat has 
the same measurements, with a negative y value. 



Patch 


Design Curve 
Control Points 


Knot 

Vector 


Endpoints 


Plane 

Definition Point 


Gaussian 

Curvature 


Forward 

Side 


(1.40, 0.00, 5.30) 
(4.17, 2.32, 4.26) 
(10.7, 5.86, 2.56) 
(18.2, 7.09, 1.87) 
(21.9, 7.38, 1.66) 


0, 0, 0, 0 
0.5 

1, 1, 1, 1 


(0.00, 0.00, 9.00) 
(20.0, 8.30, 6.74) 


(20.0, 0.00, 6.74) 


8.18 X 10-^ 


Forward 

Bottom 


(1.40, 0.00, 5.30) 
(14.0, 0.00, 0.00) 


(0.00, 0.00, 0.00) 


0.00 


Aft 

Side 


(21.9, 7.38, 1.66) 
(25.6, 7.67, 1.44) 
(33.0, 7.66, 1.28) 
(40.4, 7.46, 1.27) 
(44.1, 7.20, 1.70) 


0, 0, 0, 0 
0.5 

1, 1, 1, 1 


(20.0, 8.30, 6.74) 
(45.0, 7.70, 6.10) 


(45.0, 0.00, 6.10) 


6.62 X 10-^ 


Aft 

Bottom 


(14.0, 0.00, 0.00) 
(44.1, 0.00, 0.50) 


(0.00, 0.00, 0.00) 


9.74 X 10"'’ 



Table 8.2: Developable Surface Data, Starboard Side 
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Constrained by Plane 

All the second directrices were constrained so that they lie in a plane defined by the two 
designated endpoints of ra(ti) and a third point g provided by the user. This point g is 
termed the Plane Definition Point and is shown in Table 8.2. For the bottom surfaces, 
the constraining plane was the xz-plane which divides the boat into symmetric halves. For 
the side surfaces, the constraining plane was positioned so that a line in the plane and 
perpendicular to the x axis would be parallel to the xy plane. 

Figure 8-2 shows the resulting developable surfaces with the design lines shown as dotted 
lines for reference. Note that the forward portion of the bottom of the boat is a triangular 
degenerate patch. Figure 8-3 is a bottom view of the boat. One side shows the rulings of 
the developable surfaces while the other side is a shaded view. 






Figure 8-2: Boat Composed of Developable Degree (3-1) B-Spline Surfaces 




Figure 8-3: Boat Bottom View 



91 



Development 

The plane development of each surface of the boat was computed, thus providing cutting 
information for the flat steel. Using the method described in Section 5.4, the planar surfaces 
shown in Figure 8-4 were constructed. All developments were within a tolerance of 10“®. 




(a) Aft Side (b) Forward Side 




Geodesics 

At the Electric Boat Corporation shipyard, many developable surfaces are used in the 
construction of submarines, since a large portion of the hull is cylindrical. When bending 
the steel into shape, accuracy is checked with detailed measurements of geodesics drawn on 
the surface. When the surface is flat, straight lines are drawn between the corners of the 
surface, forming an X as shown in Figure 8-5(a). The position of this line on the curved 
surface is easily computed using the method described in Chapter 6. As an example, the 
geodesics on the Forward Side section have been calculated. The uv parameter values are 
shown in Figure 8-5(b) and the geodesics on the developable surface are shown in Figure 8- 
5(c). 

Lines of Curvature 

These surfaces have no inflection lines, but the lines of curvature were calculated as described 
in Chapter 7 and are shown in Figure 8-6. The dashed lines represent lines of minimum 
curvature and the solid lines are lines of maximum curvature. Note that one set is always 
composed of straight lines that correspond to the generators. 
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Y 





(a) Surface Developed onto Plane (b) uv Parametric Space 




Figure 8-5: Geodesics on the Forward Side Section 
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(a) Aft Side 




(d) Forward Bottom 



(c) Aft Bottom 



Figure 8-6: Lines of Curvature on Boat Surfaces (Solid Lines are Majcimum Curvature 
Lines, Dashed Lines are Minimum Curvature Lines) 
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Chapter 9 



Conclusions and Recommendations 



9.1 Conclusions 

In this thesis a novel method to design complex objects in terms of B-spline developable 
surfaces focusing on user friendliness is developed. The user only needs to specify the 
design B-spline curve and the two end rulings, and the developable surface in a B-spline 
representation is automatically generated. Although optimization techniques are relied 
upon, in most cases the computed results have a Gaussian curvature on the order of 10“^ 
or less unless the design of a very complex surface is attempted. The computational time 
for the design of all the examples in this thesis is on the order of a few tenths of seconds on 
a 100 MHz Silicon Graphics workstation. 

In addition, two differential geometry properties of developable surfaces which have 
practical engineering applications were investigated. The first property was the definition 
of inflection lines on developable surfaces, which are critical in the manufacture of such 
surfaces. We proved that an inflection line on a developable surface occurs on a generator 
which consists of entirely flat points, and derived a single variable nonlinear polynomial 
equation to locate it. It was also found that for each flat point on an inflection line of a 
developable surface, there is only one line of curvature that passes through it, and that line 
of curvature is orthogonal to the direction of the generator. 

The second differential geometry property investigated was the location of geodesics 
between two points on a developable surface. It was shown that the problem can always be 
converted from a boundary value problem to an initial value problem with a corresponding 
savings in time for the solution. 

Finally, these techniques were applied to a boat model, a complex object that was 
represented entirely with B-spline approximate developable surfaces with cubic directrices. 
The Gaussian curvature was on the order of 10"^. Geodesics and lines of curvature were 
calculated on these surfaces as well. 

9.2 Future Work 

The ultimate goal of this research is to be a useful tool in the design and manufacturing of 
ship hulls. Therefore, future research should take this work in that direction. Some specific 
areas of possible research are described below. 

Approximation. The approximation of free form surfaces by developable surfaces would 
be extremely helpful in converting existing ship designs to ones including more developable 
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surfaces. Some preliminary research has been completed by Elber and Cohen [9] into repre- 
senting surfaces as ruled surfaces, which could possibly be extended to include developable 
surfaces. In addition, research on approximation using developable surfaces was conducted 
by Hoschek and Pottmann [23] using the representation of developable surfaces as the en- 
velope of a family of planes. Alternatively, Cho [6] discussed the approximate development 
of non-developable doubly curved surfaces onto a plane, which may also be useful in the 
above approximation problem. 

Manufacturing. Once the surface is developed into a plane and the steel is cut to the 
correct size, it must be bent into the proper shape. Additional research should investigate 
the proper amount of curvature to be induced by rolling or pressing and instructions for 
the shop floor should be output. Although the lines of curvature are instrumental in this 
bending, mechanical considerations such as springback must be included in the design as 
well. Much research has been done in this area by Hardt et al [19, 18, 17]. The ultimate 
goal is for the engineer to design the ship hull and output information that is directly useful 
to the workers on the shop floor. To this end, the operating characteristics of bending 
equipment should be studied and properly accounted for. 

Geometry. Several areas of interest remain for investigation within the differential ge- 
ometry discussed in this thesis. First, a suitable non-dimensional Gaussian curvature ex- 
pression should be developed to better judge the developability of surfaces designed using 
the minimization process described in Section 5.3. 

Although strip surfaces of any degree can be handled by the computer programs devel- 
oped during this research, the programs for the surfaces with arbitrary 3D directrices can 
handle only B-spline developable surfaces with cubic directrices. It may be beneficial to 
investigate higher degree surfaces which may include more flexibility and would be able to 
more closely approximate the desired output surface. 

In addition, the surfaces in this thesis are all integral B-spline surfaces. Extension 
to include rational B-spline developable surfaces such as those investigated by Pottmann 
and Farin [37] and Lang and Roschel [29] would allow the representation of more complex 
surfaces. 

One difficulty in designing surfaces using this method is that, for surfaces with 3d 
directrices, there is no method to determine in advance if the surface will include the edge 
of regression. Therefore, regular surfaces are achieved through trial and error, and hence 
regularity conditions should be investigated more fully. 

Design. Much research hcts been conducted in design for producibility in ship hulls. For 
example. Lamb [27] shows a method to construct a bulbous bow from regular surfaces and 
discusses the simplification of bow and stern sections to include more developable surfaces. 
This should be expanded to investigate the design of an entire ship hull from developable 
surfaces. Once a ship hull is designed, it should be hydrodynamically tested to compare the 
performance to a similar hull that is not constructed of developable surfaces. The relative 
magnitude of the savings due to producibility and the losses due to performance should be 
compared. 
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